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SUMMARY 


A method for assessing the wall interference in transonic two-dimensional 
wind tunnel tests including the effects of the tunnel sidewall boundary layer 
has been developed and implemented in a computer program named TWINTN4 . The 
method involves three successive solutions of the transonic small disturbance 
potential equation to define the wind tunnel flow, the equivalent free air flow 
around the model, and the perturbation attributable to the model. Required 
input includes pressure distributions on the model and along the top and bottom 
tunnel walls which are used as boundary conditions for the wind tunnel flow. 

The wall-induced perturbation field is determined as the difference between the 
perturbation in the tunnel flow solution and the perturbation attributable to 
the model. The methodology used in the program is described and detailed 
descriptions of the computer program input and output are presented. Input and 
output for a sample case are given in an appendix. 


INTRODUCTION 

A transonic computation procedure for assessing the wall interference on an 
airfoil section tested in a two-dimensional wind tunnel was outlined in 
reference 1 and described and demonstrated in reference 2. Reference 3 is a 
user's guide to the computer program TWINTAN in which this procedure is 
implemented. Measured pressure distributions on the upper and lower tunnel 
walls and on the airfoil model are used as input and the procedure calculates 
velocity perturbations over the model due to interference of the upper and lower 
walls under the assumption that the tunnel flow is two-dimensional. 

Barnwell showed in reference 4 that the three-dimensional distortion caused 
by interaction between the pressure field around the airfoil and the thickness 
perturbations of the boundary layers on the tunnel side walls was equivalent, 
under the Prandtl-Glauert similarity rule, to a reduction in Mach number of a 
two-dimensional flow. Sewall in reference 5 reformulated the sidewall boundary 
layer effects using the Von Karman transonic similarity rule and showed that his 
formulation was capable of correlating experimental data obtained with various 
artificially thickened sidewall boundary layers. 

The study of reference 6 identified two alternative procedures for 
combining the Sewall method of correcting for the sidewall boundary layer 
effects with the TWINTAN method of assessing the interference arising from the 
upper and lower walls. In the sequential procedure, the Sewall method is used 
first to define a transformed flow free of sidewall boundary layer effects which 
is then assessed by the TWINTAN method to quantify the upper and lower wall 
effects. In the unified procedure, the existence of sidewall boundary layer 
effects is recognized in an altered governing equation used in the solution of 
the wind tunnel flow so that the assessment procedure quantifies the 
interference effects from all four walls at once. A revised version of the 
TWINTAN computer program, capable of performing the combined four-wall 



interference assessment by either of the two procedures, has been developed and 
is named TWINTN4 . The present paper provides documentation to aid the user of 
the TWINTN4 program. This paper presents the methodology used in the program 
and a detailed description of the input required and output provided by the 
program. The input data is summarized in Appendix A and input and output for a 
sample case are given in Appendix B. 

One significant difference between the two programs should be noted. In 
the TWINTAN program, the Mach number correction was determined from the local 
wall-induced velocity perturbation at a user-specified match point location. 
TWINTN4 does not use a match point. Instead, the corrected Mach number and 
angle of attack result from minimization of the mean square difference between 
the surface velocity over the airfoil in the tunnel and that over the same shape 
airfoil at the same lift in free air. In effect, this procedure is the same as 
that suggested by Murman in reference 7. In TWINTN4 , however, the minimization 
process is imbedded in and performed simultaneously with the relaxation of the 
potential solution for the free air case. Thus, the free air pressure 
distribution which best matches that in the tunnel is found with little 
additional computation time. 


SYMBOLS 

a coefficient of (f)^ in boundary condition expressions 

A,B,C,D coefficients in tridiagonal form of governing equation 
B' ,D' coefficients in line relaxation algorithm 
b half width of tunnel 

c airfoil chord 

c section lift coefficient 

l 

f transformation derivative in longitudinal direction 

g transformation derivative in transverse direction 

H nominal shape factor of sidewall boundary layer 

j longitudinal grid index, also longitudinal computational coordinate 

k transverse grid index, also transverse computational coordinate 

M Mach number 

n stretching parameter in coordinate transformation 

P,Q quantities in the finite difference form of the governing equation 

q dynamic pressure 
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s distance along airfoil surface 

5 term containing sidewall boundary layer properties 

U longitudinal velocity of uniform flow 

V resultant velocity 

x longitudinal physical coordinate 

y transverse physical coordinate 

z general physical coordinate 

a angle of attack 

P circulation 

y ratio of specific heats 

A prefix denoting an increment 

6 scaling parameter in coordinate transformation 

6* nominal displacement thickness of sidewall boundary layer 

£ general computational coordinate 

A coefficient of longitudinal term in governing equation 

p switching parameter in finite difference formulation of the governing 

equation 

$ total velocity potential 

<(> dimensionless perturbation velocity potential 

Subscripts : 

x,y denote partial differentiation with respect to x or y 

j ,k denote value at j-th or k-th grid line 

R reference condition in transformed tunnel flow 

T wind tunnel reference conditions 

co corrected far field conditions 
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METHODOLOGY 


Wall Interference Assessment Procedure 


The procedure implemented in program TWINTN4 involves three successive 
finite-difference transonic flow solutions. To allow comparison betwen 
solutions with different values of the far field velocity Uco a dimensionless 


perturbation potential <f> is defined relative to the total dimensional velocity 
potential $ by the relation 


0 


R 



x 

c 


+ 


\ 


a) 


where IL is the velocity corresponding to the reference Mach number. This 
K 

perturbation potential is governed by an extended form of the transonic small 
disturbance equation expressed as 


Atj) + <J> = 0 

xx yy 
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The tunnel sidewall boundary layer effects are introduced through the term S 
as discussed in the section headed "Sidewall Boundary Layer". The higher order 
term in ^ is included to improve the equivalence between changes in <J> x and 

Uoo- The first solution is a calculation of the flow around the test airfoil in 
the wind tunnel. For this solution U TO is set equal to U D • Pressure 

distributions measured during the test to be assessed, both on the airfoil and 
at the upper and lower walls, are imposed as boundary conditions to enhance the 
fidelity of reproduction of the actual tunnel flow. The assumptions and 
inaccuracies usually encountered in modeling the ventilated walls and the 
remaining tunnel geometry are thereby avoided. The airfoil shape defined by 
this solution is equivalent to the displacement effect of the model plus upper 
and lower surface boundary layers. 

The second solution is a calculation of the free-air flow around this 
equivalent airfoil shape. As this solution relaxes toward convergence, the 
angle of attack correction is updated to satisfy the Kutta condition for a lift 
corresponding to that in the tunnel flow, and the far field velocity and Mach 
number are updated according to an optimization algorithm which seeks to match 
the airfoil pressure distribution in free air as closely as possible to that in 
the tunnel flow. Upon convergence, the correction to angle of attack and the 
corrected Mach number are immediately available. The measured airfoil pressure 
coefficients and lift coefficient can now be adjusted to the new reference Mach 
number and the calculated pressure distribution in free air is available for 
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comparison. Differences between the two pressure distributions are attributable 
to higher order interference effects (nonuniformity of the wall-induced velocity 
field) . If these differences are judged small enough to have a negligible 
effect on the airfoil boundary layer properties, the free air pressure 
distribution can be considered to be fully corrected for wall interference. 

The third solution is another free-air solution at the corrected far field 
tfach number and is performed to define the flow perturbation attributable to the 
model. The wall-induced perturbation field is then determined as the difference 
between the model perturbation and the total perturbation defined by the tunnel 
flow solution. For this purpose, the model boundary conditions for the third 
solution are imposed in the form of singularity distributions extracted from the 
tunnel flow solution. With this model specification the airfoil camber shape is 
free to be distorted in response to removal of the wall-induced velocity field. 

Some features of the method are described in more detail in the following 
sections to aid the user in adapting the program to a specific facility as well 
as in making refinements. 


Sidewall Boundary Layer 

At the option of the user, either a sequential or a unified procedure (see 
reference 6) is used to account for the sidewall boundary layer. With either 
option, the sidewall boundary layer effects are embodied in the term S which 
appears as an increment in the coefficient of the longitudinal terra in the 
governing equation (2). This term is evaluated initially as 


/ 




+ 1 " 


“t 


The sequential procedure uses the similarity rule developed by Sewall 
(reference 5) to define an equivalent two-dimensional tunnel flow. The 
transformed reference Mach number is found by solving 


Mr ^ 

(1-M R 2 > 3/4 (1-M t 2 +S) V 4 


and all input pressure coefficients and lift and drag coefficients are then 
transformed by multiplying them by 

/ i-m t 2 +s \ ** 


x -v 

The value of S is then set to zero for all of the finite difference solutions. 
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Under the unified procedure, no similarity transformation is made; instead, 
the sidewall boundary layer effects are entered into the governing equation used 
for the tunnel flow solution through the modifying term S, and the reference 
Mach number is set equal to The value of S is then set to zero only 

for the subsequent free-air flow solutions. 



In defining the transformation derivatives f and g, the quantities j and 
k each serve the dual role of computational variable and grid index because the 
transformation is scaled so that the mesh size is one unit of the computational 
variable. The first derivatives of 4> appearing in the preceding expressions 
are all formed by single mesh differencing; for example: 
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( Vj,k+J2 g k+% (<f ’j,k+l~ <J, j,k ) 


In the following discussion, points in the lower half plane will be 
identified by a negative sign on the k index* Accordingly, g_^ = -g^- 

The governing equation is arranged in the tridiagonal form 


A A 


( 3 ) 


along j = constant lines and is solved by the method of successive line 
overrelaxation subject to boundary conditions imposed in Dirichlet form at the 
outer boundary and in various forms at the airfoil* The development of these 
boundary conditions for the three flow computations is described in subsequent 
sections. The implicit algorithm used to solve equation (3) along each line is 
based on the recursive relationship 


,k D k + B k^j , k-1 


(4) 


where the coefficients B' and D' are defined by 

n 1 — 


\ “ c k D ’k + i 
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and are calculated recursively starting with the outer boundary conditions and 
proceeding inward to k = ± 2. 

Because of the symmetric grid indexing used, the k = 1 grid line is 
treated numerically as a boundary along its entire length. At points ahead of 
the airfoil and in the wake, the boundary condition is derived from the 
governing equation. In the wake, the potential jump relation 


*3,-1 **3,l- r 


(5) 


is used along with equation (4) applied at k = ±2 to write equation (3) at 
k = 1 and obtain the solution 


*3,1 


V C 1 [ D 
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B 1+ C 1 (b ’ 2 +b '_ 2 ) 

Equations (5) and (4) are then used to complete the line solution. The same 


( 6 ) 
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procedure, but with r = 0, is used ahead of the airfoil 


Boundary Conditions for Tunnel Flow Computation 

All input pressure coefficients are interpolated onto the appropriate 
longitudinal mesh center points, converted to the form of velocity magnitude 
squared using the exact pressure coefficient definition, and stored in this 
form. The boundary conditions for the tunnel flow computation are derived 
primarily from these data. 

On the longitudinal outer boundaries at k =± KW, values of <j> are found 
from X 


V 2 + (l+<j> ) 2 +<p 2 
T x T y 


(7) 


4 > 


Values of 

y 

The <J> 


are 


and integrated to determine the boundary values of <f 
assumed initially to be zero and then are updated iteratively 

A 

gradients approaching the ends of the upper and lower boundaries are used with 
the governing equation (2) to find $ at each boundary corner. At the two 


upstream corners, 


yy 

is specified by input quantities. The variation of <f> 


on the upstream boundary is described as a fourth degree polynomial in y whose 
coefficients are determined from the values of <f> and noted above and 

y yy 

the assumption that ♦ = 0 at the lower upstream corner. On the downstream 
boundary, <J> is defined as a cubic in y with the term r/2 sgn y added to 
provide the needed potential jump across the wake. The four coefficients are 
determined to satisfy the known values of <j) and cf> at the corners. 


The airfoil boundary conditions are imposed in Dirichlet form on a slit at 
y = 0 rather than on the airfoil surface. It is not appropriate to require 
that the airfoil surface velocities be reproduced on the slit. Instead, we 
require that the longitudinal distribution of velocity potential on the airfoil 
surface be reproduced on the slit. Thus, 



(1+<|> )dx 


( 8 ) 


where V is the airfoil surface velocity determined from the pressure 
coefficients and ds is the element of surface arc length. The present program 
calculates <j> by equating the integrands of equation (8) using the 
approximation 3 ^ 


ds 

dx 


(l+<j> 

y 
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in which ♦ is calculated at the slit and is updated iteratively. This 
procedure may be recognized as an adaptation of the Riegels rule. 

The approximations inherent in the above procedure are unacceptable at a 
blunt leading edge, therefore equation (8) is used to determine <f> x for only 

those mesh intervals downstream of the leading edge. In the mesh interval 
bracketing the leading edge, upper and lower surface values of 4> x are 

determined such that two broad constraints are satisfied. The first requires 
that the total circulation determined by integrating A<f>^ over the entire 

airfoil chord must correspond to the prescribed lift coefficient. The 
circulation error is corrected directly by adjusting the Ac|> in the leading 

edge interval. The second constraint is the airfoil thickness closure condition 
which requires that the net source strength determined by integrating A(j)^ over 

the entire airfoil chord must correspond to the prescribed drag coefficient. 

This constraint is implemented by allowing the error in net source strength to 
control the average value of <|> in the leading edge interval. The direct 

effect of negative <f> x at the leading edge is to depress the value of 4> at 

the first grid point behind the leading edge, thereby creating source strength 
at this point. As a secondary effect, however, added source strength is created 
over the entire chord with a distribution such that the <f> x boundary conditions 

are still satisfied in each mesh interval. The existence of such a source 
distribution which has no effect on <|> anywhere on the chord line causes the 

thickness design problem to be ill posed without a specific constraint on 
thickness closure. The role of the thickness closure constraint in the airfoil 
thickness design problem is analogous to the role of the Kutta condition in the 
airfoil lift analysis problem. 


Solution Process for Equivalent Free-Air Flow Computation 

In the equivalent free-air computation, the equivalent airfoil shape 
determined by the tunnel flow solution is placed in a free-air environment at an 
angle of attack and Mach number such that the lift is identical with that in the 
tunnel and the distribution over the airfoil surface of total velocity magnitude 
is a best fit to that in the tunnel. The free-air environment is established by 
using outer boundary conditions calculated from a far field analytic expression 
which includes the two-dimensional lift and thickness terms given in reference 8 
and an additional term representing the far field influence of the net source 
strength at the airfoil. 

The airfoil shape is defined by values of <j> ^ evaluated at grid points on 

the slit representing the airfoil. These values may be expressed in terms of 

<f> at surrounding grid points by using a $ equal to -A<f> • Using the 

yy xx 

tridiagonal coefficients of equation (3), the results can be written 


B, <f>. 1 + C< J>. O = D 1 “ 




y 


(9a) 
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(9b) 




The values of 


evaluated by these expressions from the tunnel flow solution 


are combined with an angle of attack correction and imposed as the airfoil 

boundary condition for the equivalent free-air computation in the following form 


'j.±i 


D ±r CD '±2 1 a [ ( Vi.±i - to ] 


E ± 1 +CB '±2 


(10) 


where the upper and lower signs refer to the airfoil upper and lower surfaces, 
respectively. At each relaxation step of the potential solution, the value of 
A« is adjusted to drive the potential jump at the trailing edge toward 
agreement with the specified circulation r which is imposed as the potential 
jump across the wake. 

At certain relaxation steps, the far field velocity ratio defined as 


U 

o 

u f = u: 


is updated in accordance with the following algorithm in which the superscript 
i designates the relaxation step just completed. 

U 1+1 = U- 1 + AU 1 
f f f 
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In this process, the airfoil surface velocity magnitude V is evaluated at each 
mesh interval on the airfoil slit and the error function is formed as the 
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root mean square difference between the squares of this velocity distribution 
and that on which the airfoil boundary condition in the tunnel flow was based 
(see equation (8)). The interval Ai is the number of iteration steps since 
the last previous update and depends on the state of relaxation convergence 

but in no case is less than five. 


After each update, the following quantities are evaluated so that the 

far field boundary conditions and the potential jump across the wake may be 
updated. 
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The velocity ratio is included in the circulation update to account for the fact 
that r (like ) is expressed in units of U R c . 

The solution is considered to be converged when both a conventional 
convergence criterion for the potential relaxation is met and the magnitude of 
AEy/Ey has been below a specified level for two successive evaluations of Ey. 

After convergence, the airfoil pressure distribution in free air is calculated 
from the airfoil surface velocity magnitudes by the exact C equation and the 

experimental values of lift and pressure coefficients are corrected to the new 
reference condition by 


C 

P>°° 


= c 


£,R 


3 

V 1 1A 2 


1 n ^R 
M. 2 p,R 


% 

% 


( 13 ) 

( 14 ) 


The first term in equation (14) represents the change in free-stream static 
pressure under the assumption of constant total pressure. 

If the unified option for tunnel sidewall boundary layer was selected, 
equations (13) and (14) correct the coefficients from the tunnel test for 
interference from all four walls. Under the sequential option, however, the 
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coefficients had already been adjusted by the similarity rule for sidewall 
boundary layer and equations (13) and (14) then provide the additional 
correction for the upper and lower wall effects. 


Boundary Conditions for Model Perturbation Computation 

The purpose of the model perturbation computation is to define that part of 
the tunnel flow perturbation which is attributable directly to the airfoil. 
Accordingly , the airfoil boundary condition in this computation should assure 
that the vorticity and source distributions on the slit representing the airfoil 
duplicate those existing in the tunnel flow solution. Thus, the quantities 

A(f> . = cj> . .. - <f) . n (15a) 


(A<J> ) . = ( 4> ) . .. — ( 4> ) - i 

y j y y jj ? -i 


(15b) 


as extracted from the tunnel flow solution are to be reproduced in the free air 
computation. From equations (15), (9) and (4) we find 




Dl+D - 1 ~ C(£)t 2 +D, -2 ) + (B -1 + CB '_2 )A(j) j ~ a(A(j) y ) j 


B 1 + B -1 + c ( bT 2 + b ! _ 2 ) 


( 16 ) 


which is used along with equation (15a) to impose the airfoil boundary 
conditions in the model perturbation computation. 

In order that the calculated model perturbation represent as closely as 
possible the perturbation attributable to the model in the tunnel flow solution, 
it is important, because of the nonlinear governing equation, that the total 
Mach number environment for the calculation be similar to that in the tunnel. 

To this end, the far field Mach number used for the model perturbation 
calculation is that determined during the equivalent free-air solution so that 
the far field Mach number correction will replace, at least to first order, the 
Mach number perturbation attributable to the walls in the tunnel flow solution. 
The free-air outer boundary conditions imposed are also the same as those used 
in the equivalent free-air solution. 

After convergence, the distribution on the tunnel axis of the wall-induced 
velocities are found by subtracting the model perturbation from the total 
perturbation calculated in the tunnel flow solution. 

DESCRIPTION OF INPUT 

Program TWINTN4 is written in FORTRAN 4 and has been implemented on the CDC 
CYBER 175 computer with a central memory requirement of approximately 70,000 
(octal) 60-bit words. The program input is provided in two namelist blocks. 
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The first, labeled NPUT, contains 12 parameters to define the computational 
grid, 8 parameters for problem constraint and 8 parameters for numerical process 
control. The second block, RDPX, provides arrays of wall and model pressure 
coefficients and their longitudinal locations. All input data are 
nondimens ional with the airfoil chord assumed to be the unit of length. The 
input parameters are summarized in Appendix A and are discussed in more detail 
in the following sections. 

The program calls several subroutines from the Langley Research Center 
graphics library to create a plot vector file which can then be post-processed 
to obtain plotted output. Some modifications to the program might be required 
to obtain plots on a different computer system. 


Computational Grid Inputs 

The perturbation potential is stored by the program in an array addressed 
as PHI(J, K, L) . The index L takes on values of 1 or 2 for the upper or lower 
half plane respectively. J and K are the longitudinal and lateral indices 
respectively in each half plane. The physical configuration of the upper half 
grid is sketched in figure 1. The K = 1 grid line is the axis of grid 
symmetry between the two half planes and represents the tunnel center line, 
y = 0. The lower half grid is a mirror image of the upper. The longitudinal 
grid spacing is uniform from J1 to J4. The longitudinal spacing upstream from 
J1 and downstream from J4 and the lateral spacing of the entire grid are 
stretched by transformations of the form 


z 


n6 


C 

n-C 


(17) 


where z represents a physical coordinate and C is the computational 
coordinate in the corresponding direction. In the present program, the 
transformation is truncated to prevent z reaching infinity. Equation (17) is 
used to relate the coordinate x or y directly to the grid index J or K 
by assuming a computational mesh interval AC Q f unity. The stretching 
parameter n then represents the number of mesh intervals which would cover the 
range 0 z <_«> , and the scaling parameter 6 is the limit as z approaches 

zero of the physical mesh interval Az. The specific relations for the 
upstream, downstream and lateral transformations are formed by substituting the 
parameters tabulated below into equation (17). 


Index Range 

1 < j < ji 

J4 < J i J5 

1 < K < KB 

c 

1 

r*H 

J - J4 

K - 1 

Z 

x(Jl) - X 

x - x(J4) 

y 

n 

M 

M 

N 

6 

DXMN 

DXMN 

DYMN 
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The program sets the uniform mesh interval Ax in the range from J1 to J4 equal 
to DXMN, locates J3 at the airfoil trailing edge (x = 1), and determines J2 such 
that the airfoil leading edge (x = 0) lies in the mesh space between J2 and 
J2 + 1. The boundaries of the uniform longitudinal mesh region are 

X(J1) = 1. - DXMN * ( J3 “ Jl) 

X(J4) - 1. + DXMN * (J4 - J3) 

The foregoing information is sufficient to determine the physical 
coordinates x,y of each grid intersection J,K in terms of the following 
program input parameters. The longitudinal grid used in the free-air 
computation is defined by the integers Jl, J3, J4, J5 and M and the real 
constant DXMN which must satisfy the following constraints. 

J5 < 100 

Jl < M 

J5 - J4 < M 

( J3 - Jl) * DXMN > 1. 

J4 > J3 

The grid lines serving as the upstream and downstream boundaries of the tunnel 
flow computation are denoted by the integers JU and JD respectively which must 
satisfy: 

JU > 1 
JD < J5 

In addition, all wall boundary points, located at mesh centers from x(JU + 1/2) 
to x(JD - 1/2) must lie within the longitudinal range of the input wall 
pressures . 

The lateral grid is defined by the integers KB, N and KW and the real 
constant DYMN which must satisfy: 

KB < 50 

KB < N 

KW ^ KB 

KW identifies the lateral boundary of the tunnel flow computation and must be 
selected along with N and DYMN such that y(KW) is the lateral distance from the 
tunnel center line to the upper or lower line of input wall pressures. 

The following points should be considered also in selecting grid input 
parameters. The outer boundaries of the free-air computation, x(l), x(J5) and 
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y(KB) should lie several chord lengths from the model for best accuracy of the 
free-air boundary conditions. If, however, J1 and (J5 - J4) closely approach M, 
and KB closely approaches N, the solution convergence rate is adversely 
affected . 


Problem Constraint Inputs 

Seven real and one integer input parameters provide constraints needed in 
the tunnel flow computation. AMREF is the reference Mach number corresponding 
to the tunnel reference conditions used in reducing the input data to 
coefficient form. The lift coefficient CL defines the total airfoil 
circulation. The accuracy of the calculated wall-induced upwash depends on the 
accuracy of the CL input. CD is the airfoil drag coefficient which provides the 
airfoil thickness closure constraint. 

The parameters SLA and SLB are the onset flow streamline slopes at the 
upper and lower corners respectively of the upstream boundary of the tunnel flow 
computation x( JU) . The present formulation presumes that this upstream 
boundary is located in the solid nozzle region of the tunnel where the flow 
direction is constrained by the known wall slopes SLA and SLB. If this 
condition cannot be met or if other conditions, such as poor resolution of the 
wall pressures prevent accurate modeling of this upstream region, SLA and SLB 
may be input as zero and the tunnel flow direction monitored by comparing the 
orientation of the airfoil shape calculated in the tunnel flow solution with the 
known incidence of the test airfoil. Nominal properties of the tunnel sidewall 
boundary layer at the model location are input through DST0B2 which is the 
displacement thickness divided by the tunnel half width, and H the shape factor. 
Default values of 0. and 1. respectively are provided for these two parameters 
and define the case of no sidewall boundary layer. 

The integer parameter ISEQ, if set to 1 invokes the sequential procedure 
for accounting for the tunnel sidewall boundary layer. The default value of 
zero results in the unified procedure being followed. 


Computational Control Inputs 

The parameters discussed in this section affect the generation of a plot 
vector file, and the convergence rate, stabilization, termination and amount of 
diagnostic output for the iterative process in each of the three solutions. 
Default values of each of these parameters are established in the program but 
may be overridden simply by including new values in the NPUT namelist. 

The integer parameter IPLOT, if set to 1 causes generation of a plot vector 
file. The default value of 0 bypasses this feature. The two real parameters 
RFA and EPSU are used only in the equivalent free air solution. RFA is the 
feedback gain determining the angle of attack response to a circulation 
mismatch. EPSU specifies the convergence limit for the far field velocity 
update . 


15 


The five remaining parameters are dimensioned quantities denoting values 
stored in three-element arrays. The array subscripts from 1 to 3 correspond to 
the three successive flow solutions. The real array ORF contains the 
overrelaxation factor used at subsonic points in the successive line 
overrelaxation procedure. The real array DXT contains a damping factor to 
retard updating of longitudinal velocity perturbations. The real array EPS 
contains the convergence criterion in terms of the absolute value of the 
perturbation potential change in one iteration at the grid location where this 
change is greatest. Iteration is terminated either when this criterion is met 
or when the number of iteration steps exceeds the value obtained from the 
integer array IMX. Information describing the convergence history is available 
for output at each step of iteration. The values stored in the integer array 
ITRC determine how much of this information is actually printed. An ITRC value 
of 1 causes printing at each iteration step; a value of 2 causes printing every 
tenth step and at termination; and the default value of zero causes only the 
final values to be printed so that it may be determined whether the iteration 
converged or ran the maximum number of steps. If either of the first two 
solutions is terminated by reaching the maximum number of iterations and the 
final value of the convergence parameter is greater than 10 times the criterion 
given in EPS , then all further processing of that case is aborted. 


Pressure Distribution Inputs 

In the tunnel flow computation, the outer boundary conditions are 
established from the input pressure coefficient distributions along or near the 
upper and lower tunnel walls. Similarly, the inner boundary conditions are 
established from input pressure coefficient distribut ions on the airfoil upper 
and lower surfaces. These data are input through parameters in the namelist 
block RDPX. The wall pressure data are described in a two-element integer array 
and two real two-dimensional arrays. The integers NW(1) and NW(2) are the 
numbers of pressure input locations along the upper and lower walls 
respectively. The first real array XW inputs the x coordinates of the wall 
pressure locations in order from upstream to downstream. The first index 
numbers these locations from 1 to the appropriate NW element and the second 
index is 1 or 2 for the upper or lower wall respectively. The second real array 
CPW inputs the pressure coefficients with the same order and indexing as their 
corresponding locations. The airfoil upper and lower surface pressure 
coefficients and their chordwise locations are input in a fully ^analogous manner 
using the integers NM(1), and NM(2) and the arrays XM and GPM. All four real 
arrays are dimensioned to accommodate values of NW or NM up to 30. 

All input pressure distributions are interpolated onto boundary points 
located at the centers of the longitudinal mesh intervals. These boundary 
points range from J2 + 3/2 to J3 - 1/2 on the airfoil (K = 1) and from JU + 1/2 
to JD - 1/2 on the outer boundary (K = KW) . The interpolation routine fits a 
spline to each distribution with a zero second derivative condition at the 
spline ends. Because the routine as presently implemented does not provide for 
extrapolation, each input pressure distribution must cover the full range of 
boundary points Manual extrapolation of the experimental pressures at the 
airfoil trailing edge might be necessary in some cases to provide the required 
range of data. Near the leading edge, on the other hand, the number of 
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experimental measurements might far exceed the requirements of this program and, 
therefore, one might omit some of the experimental data. Care should be taken 
not to introduce spurious oscillations in the spline fit between data points. 

One might even input a false pressure coefficient at the leading edge to help 
the spline provide a smooth fairing in the range of interest. 


DESCRIPTION OF OUTPUT 

Output is generated during each major phase of computation so that the user 
may scrutinize the processes leading to the final wall-induced perturbation 
results. The process output will be described in four sequential phases, input 
processing, tunnel flow solution, equivalent free air solution and model 
perturbation solution. All output velocities are expressed in units of the 
reference stream velocity corresponding to the input tunnel Mach number AMREF. 
All output velocity potentials are expressed in units of the product of this 
reference velocity and the airfoil chord. 


Input Processing 

The contents of the namelist block NPUT are listed first. If the 
sequential option has been selected, the sidewall boundary layer similarity 
transformation is performed prior to all output. The transformed values of 
AMREF, CL and CD therefore appear in the NPUT listing. The longitudinal grid 
structure is next described by listing the x location and the derivative 
dj/dx. These quantities evaluated at the grid lines J are identified as X(J) 
and F1(J). The corresponding quantities evaluated at J + 1/2 are given by 
XMID(J) and F2(J). The lateral grid structure is described by the grid line 
location Y(K) , and the derivative dK/dy evaluated at the grid lines G1(K) and 
at the K + 1/2 points G2(K). The user should check to see that Y(KW) is the 
lateral location of the input wall pressures. The results of interpolating the 
airfoil pressure coefficients onto the mesh midpoints XMID(J) are then listed in 
separate groups for the upper and lower surfaces. This listing can be used in 
conjunction with the input data to judge the quality of the spline fit. 


Tunnel Flow Solution 

The convergence history of the tunnel flow solution is shown by listing the 
iteration number I, the largest change in the perturbation potential from the 
previous iteration DPHIMX, and the grid intersection JMX, KMX where this change 
is located. Negative values of KMX indicate points in the lower half plane. 

The flag ISUP is 1 if this grid point is a supersonic point and 0 if it is 
subsonic. The number of supersonic points NS and the current value of airfoil 
trailing edge thickness TTE are also listed. Both the inner and the outer 
boundary conditions depend on the solution and therefore are updated 
iteratively. To avoid using poorly developed data for the update, it is 
performed only when the convergence parameter DPHIMX is less than some threshold 
value. The threshold levels specified for the inner and outer boundary 

_3 

condition updates are 10 and EPS(1)*10 respectively. The threshold levels are 
separated so that numerical problems associated with each update process might 
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be identified more readily. Because the initial approximation to the inner 
boundary conditions is poor, it is normal for the convergence history to show 

-3 

that DPHIMX bounces up from the 10 level several times before sinking through 
it. This listing will end either with DPHIMX _< EPS(l) or with I = IMX(l) 
depending on whether or not convergence was achieved. 

The data in the following listings are identified by a heading as pertinent 
to the tunnel flow. The critical pressure coefficient is given as CPSTAR. Data 
evaluated at the airfoil surface are listed next. Two lines are printed for 
each grid index J. Data on the first line pertain to the upper surface or both 
surfaces and those on the second line to the lower surface or the increment 
between surfaces. The first four columns give data evaluated at the center of 
the mesh interval upstream from J. The calculated shape of the airfoil surface 
streamlines is given by the coordinates X(J-.5) and YMOD. The pressure 
coefficient CP(REF) can be compared with the previous listing of interpolated 
pressure coefficients to see how well the inner boundary conditions were 
satisfied. The longitudinal perturbation velocity component U is obtained 
directly from the potential gradient along the K = 1 grid line. The data in 
the last four columns are evaluated at the grid points located at X(J). V is 
the lateral component of perturbation velocity and DELTA PHI and DELTA V are the 
increments between the upper and lower surface values of the perturbation 
potential and the lateral velocity respectively. These increments are stored 
for use as model boundary conditions in the model perturbation computation. 

The perturbation velocity components along the outer boundaries of the 
tunnel flow computational region are given next. The first listing pertains to 
the upstream and downstream boundaries. The different coordinate locations 
given for the two components reflect the fact that U is evaluated at lateral 
grid lines and at midpoints in the longitudinal direction while the V component 
is evaluated at longitudinal grid lines and at midpoints in the lateral 
direction. In the case of the upper and lower boundaries listed next, the V 
component is extrapolated to grid points on the K = KW boundary. Thus, the 
listed values of U and V are given at the indicated longitudinal location along 
the upper and lower boundaries at K = KW. 

The final listing from the tunnel flow solution gives the perturbation 
potential at grid points in the immediate vicinity of the airfoil. The grid 
points are those on the K = 1, 2, and 3 grid lines in the upper and lower half 
planes for J from J1 to J4. 


Equivalent Free Air Solution 

Data output listed for the equivalent free-air solution is structured 
similarly to that described for the tunnel flow solution. This discussion, 
therefore, will highlight the differences between the two. 

In the equivalent free-air solution, the angle of attack correction DA is 
updated after each iteration for which the convergence parameter DPHIMX is less 

-3 

than 10 . The far field velocity ratio UFAC is updated either as DPHIMX 

decreases through EPS(2)*10 or if DPHIMX is less than EPS(2) provided five 
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iterations have occurred since the previous update. The convergence history 
listing includes the values of DA, UFAC and the error function EU. 

After convergence, the final far field Mach number is given as AMINF , the 
corresponding velocity ratio as UINF/UREF, the dynamic pressure ratio as QFAC, 
and the angle of attack correction is listed in radians. This far field 
condition represents the corrected reference condition for the tunnel data. 
Accordingly, the model surface pressures listed as CP( INF) as well as the 
critical pressure CPSTAR are in the form of pressure coefficients based on the 
corrected reference condition. It should be noted that the perturbation 
velocities listed in the free-air solution output are perturbations from a 
uniform flow at the far field velocity UINF. The velocity unit, however, is 
still the tunnel reference velocity or the similarity transformed velocity UREF. 

The experimental lift and pressure coefficients, corrected to the new 
reference condition, are then listed. The pressure coefficients are listed 
first at the computational mesh centers and then at the original input orifice 
locations. The root mean square increment between the calculated free air 
values and the corrected experimental values at mesh centers is given as a 
measure of the higher order interference effects not accounted for by the 
corrections to Mach number and angle of attack. 


Model Perturbation Solution 

The model perturbation solution is a simple potential relaxation with no 
boundary condition updates. The convergence history is simplified accordingly. 
After convergence, the flow properties at the airfoil surface are listed in the 
same form used for the previous solutions. 

The wall-induced velocities are calculated as the difference between the 
perturbation velocities calculated at the same locations in the tunnel flow 
solution and in this solution. This is performed in the present program only 
along the tunnel center line (K = 1) from JU to JD. The lateral velocities in 
the free air solution VF and due to the walls VW are listed at the grid point 
locations XV. The longitudinal velocity perturbations in the free air solution 
UF and due to the walls UW are given at the mesh center location XU. Over the 
model chord, values corresponding to both upper and lower surfaces are given. 
Although the wall-induced velocities should be identical on the two surfaces, 
some differences can appear in VW which reflect the accuracy with which the 
velocity jump boundary condition is satisfied. 

The final listing gives the mean of the wall-induced velocity components VW 
and UW evaluated over the model chord and the standard deviation over the model 
of these two components and their resultant. 


Plotted Output 


Two frames of plotted output are generated. In the first frame, the two 
components of wall-induced velocity are plotted against distance along the 
tunnel axis extending from about one chord length upstream of the leading edge 
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to about one chord length downstream of the trailing edge. The longitudinal 
component is plotted with + symbols and the lateral component with x 
symbols. The corrected far-field velocity ratio (relative to unity) is plotted 
with a large + at mid-chord and the angle-of-attack correction is shown by a 
large x at the trailing edge. The figure title includes the case title (first 
line of input), the standard deviations over the model of the longitudinal, 
lateral and resultant wall-induced velocities and information identifying the 
computer run, date and time. 

The second frame compares the corrected experimental airfoil pressure 
distribution, plotted with asterisks, with the calculated free-air pressure 
distribution, plotted with + symbols for the upper surface and x symbols for 
the lower surface. The most upstream points on the calculated curves are 
plotted at the center of the mesh interval bracketing the leading edge which, of 
course, might lie upstream of the actual leading-edge location. The critical 
pressure coefficient is shown by the large + symbol superimposed on the CP 
scale and the stagnation pressure coefficient is shown by the large x at the 
leading edge. The figure title includes the corrected values of Mach number and 
lift coefficient and the rms difference between the calculated and experimental 
pressure distributions. 
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APPENDIX A 


TITLE 

J 1 
J3 
J4 
J5 
M 

DXMN 

JU 

JD 

KB 

N 

DYMN 

KW 

AMREF 

CL 

CD 

SLA 

DST0B2 

H 

ISEQ 


SUMMARY OF PROGRAM INPUT 

- Title of case. Up to 80 alphanumeric characters to appear in output 

listing and plotted figure titles. 

Parameters in Namelist NPUT 

- Longitudinal grid index at beginning of uniform grid spacing. 

- Longitudinal grid index at airfoil trailing edge. 

- Longitudinal grid index at end of uniform grid spacing. 

- Grid index at downstream end of free-air domain. 

- Longitudinal grid spacing in uniform region. 

- Longitudinal grid spacing in uniform region. 

- Grid index at upstream end of tunnel flow domain. 

- Grid index at downstream end of tunnel flow domain. 

- Lateral grid indes at outer edge of free-air domain. 

- Lateral grid stretching parameter. 

- Inner limit value of lateral grid spacing. 

- Lateral grid index at outer edge of tunnel domain. 

- Tunnel reference Mach number. 

- Airfoil section lift coefficient. 

- Section drag coefficient used to define trailing edge thickness. 

- Slope of upper tunnel wall at upstream end of tunnel domain. 

- Nominal sidewall boundary layer displacement thickness divided by 

tunnel width. (Default value =0.) 

- Nominal sidewall boundary layer shape factor. (Default value = 1.) 

- Selector for sidewall boundary layer accounting method. 

= 0 for unified procedure. (Default value) 

= 1 for sequential procedure. 
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IPLOT 


RFA 

ORF(I) 

DXT(I) 

ITRC(I) 

IMX(I) 

EPS(I) 

EPSU 

NW( I) 

XW(J,I) 
CPW(J,I) 
NM( I) 

XM( J , I) 
CPM(J,I) 


Selector for plot output. 

= 0 for no plot vector file. (Default value) 

= 1 to generate plot vector file. 

Feedback gain for angle-of-attack correction. (Default value = 0.5) 

Overrelaxation factor for subsonic points in the I-th solution. 
(Default values = 1.8, 1.8, 1.8) 

Longitudinal velocity update damping factor in the I-th solution. 
(Default values = 0.5, 0.5, 0.5) 

Selector for convergence history output. 

= 0 for final values only. (Default value for all three solutions) 

= 1 for output at each iteration step. 

= 2 for output at each tenth step plus final values. 

Maximum number of iteration steps allowed in I-th solution. (Default 
values = 350, 600, 200) 

Convergence criterion for potential relaxation in I-th solution. 
(Default values = 10 10 10 ^) 

Convergence criterion for far field velocity update. (Default value = 


Parameters in Namelist RDPX 

Number of pressure input locations on I-th wall. 

1=1 for upper wall. 

1=2 for lower wall. 

Array of longitudinal locations for pressure input on I-th wall, J 
varies from 1 to NW(I). 

Array of input pressure coefficients on I-th wall, J varies from 1 to 
NW( I) • 

Number of pressure input locations on I-th airfoil surface. 

1=1 for upper surface. 

1=2 for lower surface. 

Array of longitudinal locations for pressure input on I-th airfoil 
surface, J varies from 1 to NM(I). 

Array of input pressure coefficients on I-th airfoil surface, J varies 
from 1 to NM( I) • 
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APPENDIX B 


SAMPLE CASE 

The sample case consists of the application of program TWINTN4 to results 
from a test of a NACA 0012 airfoil model in the Langley 0.3m Transonic Cryogenic 
Tunnel with slotted top and bottom walls. The ratio of tunnel height to airfoil 

chord was 4. The airfoil was set at 2.025° incidence and tested at a reference 

Mach number of 0.701 and a Reynolds number of 3.05 x 10^. 

The sequential option for accounting for the sidewall boundary layer was 
chosen for the sample case. The input card images are listed and followed by 
the program output listing. The plotted output is shown in figures 2 and 3. 
Figures 4 and 5 show the plotted output from the same case but using the unified 
procedure for accounting for the sidewall boundary layer. 
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INPUT F D P SAMPLE CASF 


NACA 0012, M= • 701 , ALPHA=2.025, R N = 3 . 0 A 6 E + 0 6 • 

SN PUT Jl r 25,J3=47,J4=49,J5=73,M=35,DXMN=.0499, 
JU=1,JD=72,KB=22,N=30,DYMN= # 05098,KW=18, 
S.LA = 0.,SLB = 0.,DRF(1)=2*1.80,ITRC<1)=2,1,2, 


IMX (1 ) = 3 50# 500,200, IPL0T=1,ISEQ=1, 

A M R P F = .701, CL = .2204,00= . 00 76 , H = 1 . 5 0 A 2, D 5 T OR 2 
*RDPX N M ( 1 ) = 24, 2 4 , N y ( 1 ) = 26, 28, 


XW(1, 1) = -.40567F+01, 

- • 3 72 33E + 01 , — *339 00 E +01, 
-•17233E+01, - . 1 3900F + 0 1, 
•27667E+00, .61000F+00, 

• 22767E + 01, • 26 1 00 E +0 1 , 

•4 276 7F +01, 

CPW{ 1,1)= - • 5 4 86 1 E- 0 2 , 
-.12283E-02, -.98122E-03, 

• 1 3 5 1 6 F-02 , - • 371536-02, 
-.60229E-01, -.59275E-01, 
-.280566-01, -.32320E-01, 


- • 30 567E+01, 
- #10567 E +01, 
•94333F+00, 
• 2 94 33 E + 0 1 , 


• 6 2 4 3 1 E-0 3 , 
-.04227E-O2, 
-.50005E-01, 
-.36963E-01, 


27233E+01, 
72333E+00, 
. 12767E+01, 
•32767E+01, 


-. 52061E-03, 
19627E-01, 
-.50088E-01, 
-•37117E-01, 


-.39552F-01, 
XW(1, 2 ) = - 
-.37233E+01, 
-.13900E+01, 
.61000E+00, 

• 267 33F +01, 

• 45067E + 01, 
CPW( 1,2 )= - 

• 691 73E-03 , 

• 2 2097E-01 , 

. 20283F-03, 

— •92554 E -02 , 

-.18373E-01, 

XM(1, 1)= . 

. 14361E-01, 
•20155F+00, 
•50149E+00, 
•80091F+00, 
CPM( 1, 1) = 

32689E+00, 

- • 8 47 2 6E + 00 , 
-.37318E+00, 

- • 740 6 2 E-01 , 

X M ( 1 , 2 ) = . 

• 11718 E — 01, 
•19846E+00, 
•49R06E+00, 
•79R59E+00, 
CPM( 1, 2 ) = 
•49190F+00, 

-•38417E+00, 

23905F+00, 

— • 4 08 1 4E — 01 , 
SEND 


.405676+01, 

-. 33900E+01, 
10567F+01, 
•943336+00, 
. 30067 E+01, 

• 48 4 00 E +0 1 , 
. 2 2 372 E-0 1, 
-•35306E-02, 

• 18 1 2 3 E— 0 1 , 
. 1417 5 E -02 , 

-.15335E-01, 
-. 17930E-01, 
2 2 477F-0 3, 

• 2 8 2 2 1 F — 0 1 , 

• 25228F +00, 
•55082E+00, 
. 85039F+00, 

. 10741F+01, 
-.61743F+00, 
-•70543E+00, 
-. 31092E+00, 
-• 20218 F — 0 1 , 
2 2477F-03, 

• 2 4 3 4 1 E -0 1 , 
•24811E+00, 
. 5 4 8 6 9 E +00, 
. 84924E+00, 

. 10741F+01, 

. 17864 F +00, 
-.37348E+00, 
-•20599F + 00, 
-.R3148E-02, 


- • 30 5 67 E +0 1 , 
-•72333 E + 00, 

• 1 27 67 E + 0 1 , 
•334 00 E+01, 
•51733F+01, 

• 47236E-02, 

• 2 2969 E — 01 , 
•84934 E-03, 

-.16095 E-0 1 , 
— • 4 5 1 4 2 E-0 1 , 

• 5 2 9 7 0 E-0 1 , 

• 30200F + 00, 
•60089E+00, 
•89945E+00, 

-.83530E+00, 

-.60428F+00, 

-•26626E+00, 

•44718E-01, 

• 4 99 36 E-0 1 , 
•29807F+00, 

• 59798E + 00, 
. 89914E+00, 

-• 10019^+00, 
-.34969 E + 00, 
-.17514E+00, 

• 29585 E-0 1 , 


-•27233E+01, 
-. 3 9000 F + 00, 
. 16100E+01, 
. 3673 3E +01, 


. 14234F-01, 
•97888E-02, 
-• 88307E— 03, 
— .18531E-01, 


• 7 797 7E-0 1 , 
. 35197E+00, 

• 65075E + 00, 
. 94878E+00, 

-•92300E+00, 
5385 7E+00, 
-•22483E+00, 
. 1 204 1 E + 00, 

• 74166F-01, 
. 34861E+00, 
•64853E+00, 
. 94766E+00, 

-•21237 E +00, 
32701E+00, 
-# 15402E+00, 

• 11052F +00, 


= 1 • 543E-02, $ 


- • 2 3900E +0 1 , 
-• 39000 E+00, 
•16100E+01, 
•361 00 E+01, 


•86244E-03, 

— • 3698 9 E—0 1 , 

— • 32 1 47 E-0 1 , 
- • 36747E-0 1 , 


— ,205 67 E+01, 

— • 5 6667E-0 1 , 

• 20067E +0 1 , 

• 384 OOF +01, 


• 98306E-02, 

— *56598 E-02, 

— .23261F-0 2, 

— • 1 8 062 E—0 1 , 


•10296E+00, 
•40145F+00, 
• 70067 E + 0 0, 
.987503 , 

- • 97980E +00, 
-.486 97 E+00, 
-.176026+00, 
.21364 , 

•98549F-01, 
•39795F+00, 
•69885E+00, 
.987503 , 

-.294136+00, 
-•30194E+00, 
11759E+00, 
.255044 , 


•20567E+01, 
. 5 66 67 E-0 1 , 
•19433F+01, 
•39433E+01, 


. 11078E-03, 
•48082F-01, 
• 332 0 4F-01 , 
•39867F-01, 


. 17233F+01, 
• 27667E+00, 
•23400E+01, 
•41733E+01, 


• 13 4 2 4 E-0 1 , 
•21510E-02, 

• 7041 8 E-02, 

• 18691 E-0 1 , 


. 15231E+00, 
•45146E+00, 
• 74982E+00, 


• 94881E+00, 
•42540E+00, 
. 12971E+00, 


. 14860F+00, 
•44826E+00, 
• 74836E+00, 


. 36746E+00, 
. 26462F+00, 
• 8576 3 E-0 1 , 
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INPUT 


SAMPLE CASE OUTPUT 


J1 

X 

25, 

J 3 

X 

47, 

J4 

• 

4 9, 

J 5 

X 

73, 

M 

X 

35, 

DXMN 

X 

•499E-01, 

JU 

3 

1, 

JD 

X 

72, 

KB 

X 

22, 

N 

X 

30, 

0 YMN 

X 

•5098E-01, 

KU 

3 

18, 

AMREF 

= 

•68711860072509E+00, 

CL 

X 

• 22335850000608E+00, 

CD 

3 

.77020172415869E-02, 

SLA 

X 

o 

• 

o 

SLB 

X 

•v 

o 

• 

o 

RF A 

3 

• 5E+00, 

ORF 

= 

• 18E + 01, • 10E + O1, .18E+01, 

DXT 

X 

• 5 E + 00, • 5E + 00, * 5 E +00, 

ITRC 

3 

2, 1, 2, 

FPS 

S 

• IE-05, .IE-04, • 1 E — 0 5 , 

I M X 

= 

350, 500, 200, 

FPSU 

s 

• IE-03, 

ISFQ 

= 

1, 

DSTOB2 

X 

.1543F-01, 

M 

X 

• 1 5042 E + 01 , 

i Pirn 

3 

1, 



J 

«ACA 001 2 » M = • 7 C 1 » 
X(J) 

ALPHA=2.025, 
FI ( J ) 

P^*3.046E+06. 
F2C J ) 

X MID ( J ) 

1 

-. 390P36+C1 

•19795E+01 

.216356+01 

-.367726+01 

2 

-.344536+01 

.23557F+01 

.255616+01 

-.324976+01 

3 

-.305346+01 

•27647E+01 

• 2961 5 E + 0 1 

— .28857E+01 

4 

-.27176F+G1 

•32064E+01 

•34395E+01 

-.25722E+01 

5 

-.242656+01 

•36808E+01 

•39303E+01 

— . 22 993 E +0 1 

6 

-.217ieF+01 

.41B80E+01 

.44538E+01 

-.20595E+01 

7 

-.1Q470F+01 

.47278E+01 

• 50 100 E +01 

-.18472E+01 

e 

-. 17473F+01 

• 5 3 004 E + 0 1 

• 55 990E + 01 

-• 1650OE+O1 

9 

-.156P5E+01 

.59057E+01 

• 6220 6E+0 1 

— • 14882E + 01 

10 

-.14G77F+01 

.65437F+01 

•68750E+01 

13349E+01 

u 

— .12621E+C1 

•72144E+01 

.75621E+01 

-.119606+01 

12 

-.1129PF+01 

. 79 1 7 9 E +0 1 

.82619E+01 

-. 10694E + 01 

13 

-.1009CF+C1 

.P6540E+01 

• 90344E+01 

-.953676+00 

1* 

-. 69026 E +00 

.942296+01 

•98196E+01 

-.847366+00 

15 

-. 7964CE+00 

•10225E+02 

, 10638F+02 

-. 74940E+00 

16 

-. 7 0 2 26 F +00 

•11059E+02 

•11488E+02 

-.658846+00 

17 

-.61 526E+0C 

•11926E+02 

. 12372E+02 

-.574876+00 

IB 

-. 53443E+C0 

•12826E+02 

. 13288E + 02 

496806+00 

19 

-.459146+00 

.13758E+02 

. 1 42 3 7 E + 02 

-•42402E+00 

20 

-. 3RBPRE + 00 

.147236+02 

. 15218F+02 

-.35603E+00 

21 

-. 32315F + 00 

.157216+02 

.162326+02 

-•29235E+00 

22 

-.2O153F+00 

. 16752E + 02 

•17279E+02 

-.232 60 E +00 

23 

-.20365F+CC 

.17*156+02 

.18359E+02 

-. 17641E+00 

24 

149176 + 00 

.16911E+02 

.194726+02 

-.123496+00 

25 

-.970OCE-O1 

.20040E+02 

•20040E+02 

-•72850E-01 

26 

-.47900F-01 

.20040E+02 

.200406+02 

— .2295 0 E-0 1 

27 

.2JCCC6-C2 

.2 0040E+02 

• 2 004 0 E + 02 

• 2 6 9 50 E-01 

2 P 

. p 1 Q 00 E -0 1 

•20 G4 06+02 

.200406+02 

.768506-01 

2Q 

. 101006+00 

.20040E+02 

.200406+02 

• 1 2 67 5E +00 

3G 

•15170F+0C 

.2C040E+02 

•20040E+02 

. 1 7 665 E + 00 

31 

.20160E+CC 

.20040E+02 

•20040E+02 

•22655E+00 

32 

. 25 1 50F +00 

•20C40E+02 

. 20040E+02 

•27645E+00 

33 

. 30 1 40 E +00 

.200406+02 

•20040E+02 

•32635E+00 

34 

. 351 3C+'+C0 

• 2004 OF + O 2 

.200406+0 2 

. 3 762 5E +00 

35 

. 4Q120E+0C 

• 20040E + 02 

. 20040E+02 

•42615E+00 

36 

.451 1CE+C0 

•20040E+02 

.200406+02 

.476056+00 

37 

. 501CCF+00 

• 20040E + 02 

• 20040E+02 

. 52595E+00 

33 

. 55C90F+G0 

.200406+02 

• 2004 0 E + 02 

.575856+00 

39 

. cOG 60 F +00 

.200406+02 

.200406+02 

.625756+00 

40 

.65C70E+CC 

. 20040E +0 2 

• 2004 0E + Q2 

• 67565E +00 

41 

. 7006C6+CC 

.20C40E+02 

.20040E+02 

•72555F+00 

4 ? 

. 7505CE+00 

.200406+02 

• 2 004 0 E + 0 2 

.775456+00 

43 

.30040F+C0 

.200406+02 

• 20040E+02 

• 82 535 E + QO 

44 

. 0 5 030 F +00 

.200406+02 

• 2004 0 E + 02 

•87525E+00 

45 

.9002CF+00 

• 20 040E + 02 

•20040E+02 

• 92 5 1 5 E + 00 

46 

• 95C1CE+GC 1 

.200406+02 

.200406+02 

•97505E+00 

47 

. 100CCF+C1 

.200406+02 

.200406+02 

• 10250E+01 

4« 

.10499E+01 

.200406+02 

. 2 004 0 E +02 

.107496+01 

40 

. 10990 6 + 01 

•20040F+02 

•19472E+02 

. 11255E+01 

50 

.1L512F+01 

.189HE + 02 

.183596+02 

.117646+01 

51 

.120666+01 

♦17815F+02 

. 17279E+02 

.123466+01 

52 

. 126356+01 

. 16 7 5 2 E + 0 2 

. 1&232E+02 

. 12943E+01 

53 

. 13252E+C1 

•15721E+02 

• 1521BE+02 

. 135P0E+01 

54 

. 1 39C9E+C1 

.14723E+02 

.142376+02 

. 142606+01 

55 

. 14611F+01 

• 13 759 E + 0 2 

. 13288E+02 

.149886+01 

56 

.153646+01 

.12B26E+02 

. 12372F+02 

.157686+01 

57 

.161736+01 

.119266+0 2 

. 11488F+02 

. 1660PE+01 

58 

. 17044E+01 

•11059E+02 

. 10638E+02 

• 17514E+01 

59 

. 17QP4E+01 

•10225E+02 

. 98196E+01 

.184936+01 

60 

. 190C3E+01 

.942296+01 

.903446+01 

. 1955b F. + 01 

61 

•20110F+01 

•86540F+01 

•B2819F+01 

.207146+01 

62 

•2131PF+C1 

.791796+01 

.756218+01 

. 21979E+01 

63 

• 2264 1 E + 01 

.72144E+01 

.68750E+01 

•23369E+01 

64 

.240976+01 

.65437E+01 

•62206E+01 

.249016+01 

65 

•257C5F+G1 

.59057E+01 

. 55990E+01 

.265986+01 

66 

. 274936+01 

. 53004 E +0 1 

•50100E+01 

• 2 84 91 E +0 1 

67 

. 2949CE+G1 

• 4727 8E +0 1 

.445396+01 

•30613E+01 

68 

. 3173PF+01 

•41P80F+01 

. 39303E+01 

. 33010E+01 

69 

. 342*56+01 

. 36808E+01 

•34395E+01 

• 3 5738 E + 01 

70 

. 37196F+01 

.320646+01 

.299156+01 

.38873E+01 

71 

. 40554F+01 

.276476+01 

•25561E+01 

.425106+01 

72 

• 444 7 3 E +0 1 

.235576+01 

•21635E+01 

•467P4E+01 

73 

.491C3F+01 

.197956+01 

• 180 36 E + 01 

.510766+01 


26 



K 

Y<K ) 

G1(K) 

G2 t K ) 

1 

0. 

• 1961 6E + 0 2 

* 1 896 7E + 02 

2 

• 527 38 E-01 

• 18 3 30E + 02 

•17703E+02 

3 

• 1 0924 F +00 

.17087F+02 

•16482F+02 

4 

.16993E+00 

• 15889F + 02 

.15306F+02 

5 

• 23529F+00 

•14733F+02 

• 1 417 2E + 02 

6 

• 3058PE +00 

• 13622F +02 

. 1308 2E + 02 

7 

• 3823 5E +00 

•12554F+02 

• 12036E+02 

8 

.46547F+00 

• 11 530E + 02 

.11034E+02 

9 

• 5561 5F +00 

• 10549F + 02 

• 1 007 5 E + 02 

10 

.65546E+00 

.96116E+01 

•91594E+01 

11 

•76470E+00 

• 07 1 80E + 01 

• 82876E+01 

12 

* 88544F+00 

* 78680E +01 

• 7459 4E + 01 

13 

.10196E+01 

• 70 61 6E + 0 1 

• 6674 7E + 01 

14 

.11695E+01 

•62908E +01 

• 59337E + 01 

15 

•133P2F+01 

• 55795E+01 

• 5236 3E + 01 

16 

.15294F+01 

• 49039F. +0 1 

• 4562 4E + 01 

17 

• 1 74 79 E +01 

•42718E+01 

• 3972 1 E + 01 

18 

• 20000F +01 

• 36 0 34 E +0 1 

• 3 405 5 E + 01 

19 

• 22941 F +01 

• 31 385E +01 

• 2082 4F + 01 

20 

• 26417E+01 

• 26372F+01 

• 24029E + 0 1 

21 

•30588E+01 

.21795F+01 

• 1967 OE + 01 

22 

• 35 686E +0 1 

•17654 E+01 

• 15747E+01 


AIRFOIL 

CP INTERPOLATED 

AT MESH 

MIDPOINTS 







-•62005 

-.93310 

-.99806 

-.91552 

-.78805 

-.65884 

-.57667 

-.52058 

-.46253 

-•<>0533 

-•34669 

-.29024 

-.24901 

-.20338 

-•15500 

-.10357 

-.04889 

.01203 

.08102 

.16388 

•13984 

-.22510 

-•35281 

-.38578 

-•38582 

-.36527 

-.34118 

-.31871 

-.28447 

-.25318 

-•22485 

-•18997 

-.16652 

-.13684 

-.10243 

-.06243 

-.02272 

.00925 

.06270 

.20789 

I JMX 

KMX ISUP 

NS 

DPHIMX 

TTE 







10 

51 

-1 

0 

0 

•30275E-01 

0 . 

20 

54 

-1 

0 

0 

• 3 7793F-02 

0 . 

30 

47 

2 

0 

1 

- . 3 192 7E-02 

-•37203F-01 

40 

25 

1 

0 

0 

- . 88 577E -03 

- • 2 1 847E-01 

50 

27 

-1 

0 

0 

-♦ 57920E-03 

• 3 3 738E-02 

60 

27 

1 

0 

0 

- • 2 2779E-03 

• 2 9849E-02 

70 

20 

-3 

0 

0 

-.12368E-03 

• 41 995E-02 

80 

22 

-5 

0 

0 

-• 57585E -04 

.41 379E-02 

90 

17 

-5 

0 

0 

-.24974E-04 

• 39589F-Q2 

100 

18 

-5 

0 

0 

-• 11480E-04 

• 38857E-02 

no 

71 

17 

0 

0 

.1674 IE-04 

. 38 654E-02 

120 

55 

-5 

0 

0 

-.43619E-05 

• 3 6 502 F-02 

130 

46 

-9 

0 

0 

-.22216E-05 

• 3 8 51 1 E-02 

140 

36 

-11 

0 

0 

-.11129E-05 

• 3 fl 528E-02 

142 

34 

-11 

0 

0 

-.96593E-06 

• 38520E-02 


ro 



RESULTS IN TUNNEL FLOw 


STAR = -.P3370E+00 


XU-.5) 

Y M 00 

CP(kEF) 

U 

-•02295 

0.00000 

.89501 

-.53092 


0.00000 

.80344 

- • 4 7 57P 

.02695 

.03972 

-.62026 

.46060 


“.02403 

.13740 

-.00508 

.07665 

.05243 

-.93314 

.46314 


-.03796 

-.22524 

.13880 

.12675 

.05961 

-.99694 

.47336 


-.04702 

-.35289 

.18311 

.17665 

.06326 

-.91519 

.42854 


“.05300 

-.38578 

.19102 

.22655 

.06458 

-.78800 

.36919 


-.05691 

-.38581 

.18755 

• 2 7 64 5 

.06393 

-.65896 

.31101 


-.05924 

-.36526 

.17640 

.32635 

.06215 

-.57676 

.27436 


-.06044 

-.34117 

.16455 

.37625 

.05959 

-.52065 

.24968 


-.06072 

-.31871 

.15390 

. A 2 61 5 

.05623 

-.46261 

.22416 


-.06014 

-.28448 

.13798 

. 47605 

.05216 

-.40540 

.19889 


-.05893 

-.25319 

.12349 

.52595 

.04751 

-.34675 

.17265 


-.05720 

-.22486 

.11043 

.57585 

.04237 

-.29027 

.14687 


-.05495 

-.18998 

.09414 

.62575 

.03697 

-.24983 

.12846 


-•05241 

-.16652 

.08327 

.67565 

.03116 

-.20340 

.10718 


-.04944 

-.13684 

.06945 

.72555 

.02500 

-.15501 

.08467 


-.04606 

-.10243 

.05325 

.77545 

.01052 

-.10358 

.06032 


-.04231 

-.06243 

.03406 

.02535 

.01177 

-.04890 

.03392 


-.03832 

-.02272 

.01463 

.87525 

.00481 

.01203 

.00376 


-.03423 

.00925 

-.00091 

.92515 

-.00226 

.08101 

-.03117 


-.02969 

.06269 

-.02655 

.97505 

-.00949 

.18384 

-.08589 

1.02495 

-.02400 

-.01601 

-.01997 

.20786 

-.10300 


J 

X< J) 

V 

DELTA PHI 

DELTA V 

27 

.00200 

.79599 





-.48157 

-.00275 

1.27756 

28 

.05190 

.25477 





-.27912 

.02049 

.53389 

29 

.10180 

.14388 





-.18164 

.03667 

.32552 

30 

.15170 

.07316 





-.11979 

.05115 

.19295 

31 

.20160 

.02642 


‘ 



-.07831 

.06301 

.10472 

32 

.25150 

-.01312 





-.04682 

.07207 

.03370 

33 

.30140 

-.03554 





-.02388 

.07879 

-.01166 

34 

.35130 

-.05135 





-.00578 

.08427 

-.04557 

35 

.40120 

-.06733 





.01172 

.08905 

-.07905 

36 

.45110 

-.08126 





.02424 

.09335 

-.10550 

37 

.50100 

-.09354 





.03458 

.09711 

-.12812 

38 

.55090 

-.10304 





.04508 

.10021 

-.14812 

39 

.60080 

-.10820 





.05094 

.10284 

-.15914 

40 

.65070 

-.11632 





.05960 

.10510 

-.17592 

41 

.70060 

-.12347 





.06777 

.10698 

-.19123 

42 

.75050 

-.12992 





.07516 

.10655 

-.20507 

43 

.80040 

-.13523 





.07987 

.10986 

-.21509 

44 

.85030 

-.13955 





.08200 

.11082 

-.22155 

45 

.90020 

-.14169 





.09101 

.11106 

-.23271 

46 

.95010 

-.14482 





.11387 

.11083 

-.25869 

47 

1.00000 

-.13079 





.08092 

.11168 

-.21171 



tJUTPP BOUNDARY VELOCITIES 



Y = 

-.36772E+01 

YU 

Y V 

U 

0. 

• 2 636 IE -01 

-.16305E-02 

. 52738F-01 

• 60 9 82 E -0 1 

-.15756E-02 

• 10Q24E+00 

. 13958F +00 

-•15 158E-0 2 

• 16993E+00 

.2 0260E+00 

-.14497F-02 

• 2 3 5 2 9E + 00 

•27057E+00 

-•13 770E-0 2 

. 30588F+00 

• 3441 OE +00 

- • 1 2 97 1 E-0 2 

• 38235F+00 

•42389P+00 

-•12092E-02 

• 465 47 F + 00 

. 51079E+00 

-.1112 4 E-0 2 

• 5 56 1 5F +00 

•60577E+00 

10061E-02 

• 6 5 5 4 6F + CO 

•71005E+00 

-.88930E-03 

• 764 70E + 00 

• 82 50 3E + 00 

- • 76 1 1 7E-0 3 

• 8 8 5 44E + 00 

•95247E+00 

- • 62 08 7F-0 3 

• 10196E+01 

•10945F+01 

-.46761E-03 

•11695E+01 

•12538E+01 

- • 30067F-0 3 

• 13382E+01 

. 14337E+01 

- • 1 1 9 3 5 E-0 3 

.15294E+01 

. 16385E+01 

• 77 3 54E-0 4 

. 17479E+01 

. 18738E+01 

•29184E-03 

0. 

-.26361F-01 

- • 1 6 30 5E-02 

-. 52738F-01 

-. 80982E-01 

16837E-02 

-• 1 092 4 E +00 

- • 1 39 58E +00 

- • 1 7 38 7E-0 2 

-. 16993E+00 

- • 202 60E +00 

-.17 95 4E-02 

-• 2 3 529E + 00 

- * 270 5 7E +00 

-. 185 3 3E-02 

-. 30586E+00 

-.34410E+00 

-•19 1 17E-02 

-.3H235E+00 

-.42389E+00 

- • 1 9696 E-0 2 

-• 4 6 5 4 7F + 00 

- • 51 079E +00 

- • 202 5 4E-02 

55615E+00 

- • 60 577E +00 

-. 20 768 E-02 

-.65546E+00 

- • 7 1 00 5E +00 

-.21 203 E-0 2 

-• 76470E +00 

-.R2503E+00 

-.21 50 8 E-02 

-. 88544E+00 

- • 952 47F +00 

-.21606E-02 

- • 1 0196E +0 1 

10945E+01 

-.21 30 0E-02 

-• 1169 5E +01 

-.12538E+01 

- • 2 0656E-02 

-. 13382E+01 

14337E+01 

- • 19173E-02 

15294E+01 

-.16385E+01 

- * 1 6 55 3E-02 

-. 17479E+01 

-. 18738E+01 

-.12 31 OE-O 2 


39083E+01 

.42510E+01 

. 44473E+01 

V 

U 

V 

-. 11308E-02 

• 1 2835E-01 

-.46250E-03 

- • 11 107 F-0 2 

• 1 2969E-01 

- • 50906E-03 

-.10874E-02 

. 1 3117E-0 1 

-.53701E-03 

10605E-02 

. 1 3279E-01 

-.56645E-03 

-. 10293E-02 

• 1 3458E-01 

- . 59747E-03 

-. 993 38 E-03 

• 1 3658E-01 

-.63018E-03 

-.95196E-03 

• 1 3880E-01 

-.66469E-03 

- • 90430F -03 

• 14129E-01 

-.70110E-03 

-.84951 E-03 

. 14409E-01 

-.73948E-03 

78666F-03 

.14726E-01 

-.77994E-03 

-. 7147 2 E-0 3 

. 15088E-01 

-.82250E-03 

- * 63269 F-03 

• 1 5502E-01 

- . 86715 E-03 

-.53 96 4E-03 

• 1 5981 E-01 

-.91381E-03 

-•43502E-03 

• 16537E-01 

-.96224E-03 

-. 31900E-03 

• 1 7191 E-01 

-. 10120E-02 

-. 19335E-03 

.17964E-01 

-. 10622E-02 

-.62918E-04 

• 18888E-01 

-. 11114E-02 

-. 11486E-02 

• 12835E-01 

- .45641E-03 

-. 11654E-02 

• 12704E-01 

-.42892E-03 

-. 11813E-02 

. 12566E-01 

-.39888E-03 

11960E-02 

• 12423E-01 

- • 36594E-03 

-. 12089E-02 

.12273E-01 

-.32969E-03 

-. 12191E-02 

• 12117F-01 

-.28963E-03 

-.12256E-02 

•11954E-01 

-.24515E-03 

-.12270E-0? 

• 11784E-01 

- • 19 553E-03 

-. 12213E-02 

• 11607E-01 

-.13987E-03 

- • 1 2059 E-02 

.11422E-01 

-.77076E-04 

-• 1 1 772 E-02 

.11231 E-01 

-.57597E-05 

-.11 300E-02 

• 11033E-01 

. 75 815E-04 

-. 10572E-02 

. 10827E-01 

. 1 6987E-03 

- • 948 28 E-03 

• 10614E-01 

. 27930E-03 

-. 78023 E-03 

• 10393E-01 

• 40791 E-03 

- • 55447 E-03 

• 10159E-01 

• 56080E-03 

- • 21 269 F-03 

. 99053E-02 

•74496E-03 



XU 

• 367726+01 
.324976+01 
. 268576+01 
• 25722E+C1 
• 229936+01 
•205956+01 
.1R472E+01 
. 165F0E+01 
.148626+01 
.133496+01 
. 119606+01 
.106946+01 
.95367P+00 
. 84736E+00 
. 749406+00 
.658846+00 
•574*76+00 
.496606+00 
. 42402E+00 
•35603E+CC 
•29235 c +00 
.232606+00 
. 176416+00 
.1234QF+00 
.728506-01 
.229506-01 
.269506-01 
. 76H5CF-01 
.126756+00 
. 176656 + 00 
. 22655E+CC 
.276456+00 
.326356+00 
.376256+00 
.426156+00 
.476056+00 
.525956+00 
. 575*56 + 00 
. 62575E + 0C 
.675656+00 
. 72555E + 00 
. 775456 + 00 
.825356+00 
.875256+00 
.925156+00 
.975056+00 
. 1025CE + C1 
. 10749E + 01 
. 112556 + 01 
. 117846 + 01 
•12346 c +01 
. 129436 + 01 
.135806+01 
. 1426CE + 01 
. 149686 + 01 
.157686+01 
. 16608 c +01 
.175146+01 
.1R493E+01 
. 19556E + 01 
.207146+01 
.219796+01 
•23369E+01 
.249016+01 
•26598E+G1 
.284916+01 
.306136+01 
•33010E+01 
•35738E+01 
, 38873E+01 
.425106+01 


Y * 

XV 

— . 39063E+01 
-.34453E+01 ‘ 
-.305346+01 
-.27176E+01 
-.242656+01 
-.217186+01 
-. 19470E+01 
-. 17473E+01 
— . 15685E+01 
-.140776+01 
-.126216+01 
-.112986+01 
-.100906+01 
-.098286+00 
-.79640E+00 
-.702366+00 
-.615286+00 
-.534436+00 
-.45914E+00 
38888E +00 
-•32315E+00 
-.261536+00 
-. 20365E+00 
-.14917E+00 
-.978006-01 
-.479006-01 
• 20000F— 02 
.519006-01 
. 101806+0C 
. 15170E+00 
.20160E+00 
.251506+00 
.301406+00 
.351306+00 
.401206+00 
.451106+00 
.501006+00 
.550906+00 
.600806+00 
.650706+00 
.700606+00 
.750506+00 

• 80040E + 00 
.650306+00 
. 900206+00 
.950106+00 
.100006+01 
.10499E+01 
. 109906+01 
.115126+01 
. 120566+01 
.126356+01 
.132526+01 
. 13909E+01 
. 14611E+01 
•15364E+01 
.161736+01 
. 17044E+01 

• 179P4E +01 
.190036+01 
.201106+01 
.213186+01 
,226416+01 
.24097E+01 
.257056+01 
.274936+01 
.294906+01 
.317386+01 
.342856+01 
.37196E+01 
,405546+01 


.20000E+01 

U 

•53134E-03 

.124576-03 

-.15474E-04 

— • 29 54 2 E-04 

— .3685 9 E— 0 3 
-.550586-04 

— .65282 E— 0 3 
- • 39 98 2 E — 0 3 

• 1 0077E-0 2 

• 22 84 1 E — 0 2 
.31613E-02 

• 41378 E— 0 2 
.54996E-02 
.724246-02 
. 9 2 966 E— 0 2 
.115056-01 
.139036-01 

• 1604 5E—0 1 

• 1 7 8 8 3 E-0 1 
.19348E-01 
.204936-01 

• 2 1 44 2 E — 0 1 

• 22 28 8 E-0 1 
.230966-01 

• 23 92 7 E — 0 1 

• 24 8 3 7E— 0 1 
.258236-01 
,268406-01 
.278396-01 
.287746-01 
.295976-01 
.30 26 3 E— 0 1 
.307326-01 
.310026-01 
.31077E-01 
.309636-01 
.306646-01 
.301866-01 
.295346-01 
.287386-01 
. 27869E-0 1 

• 27 002 E-0 1 
.262 1 1E-0 1 

• 25 572E-01 
,251596-01 
.250416-01 
.251726-01 

• 2 5 4 2 7E— 0 1 
.256846-01 

• 2 5 806 E-0 1 

• 2 5606E-0 1 

• 24 82 4E-0 1 
« 23 297E-0 1 
.211986-01 
.18852E-01 

• 1 6759E-0 1 
. 15622E-01 
.156336-01 

• 16 2 0 OE-01 
.166746-01 
.159486-01 
.146356-01 
•14058E-01 
. 1502 3 E-0 1 
.167426-01 

• 18 1 7 8 E-0 1 
. 18846E-01 

• 1 6 6 5 4 E-0 1 
.184386-01 

• 1 9 876 E -0 1 
.19993E-01 


.200006+01 

V 

.44279E-03 

• 762 1 6 E— 03 
.113016-02 

• 1661 6E-02 

• 16842E-02 
.233356-02 
.188446-02 

• 11 08 6 E— 02 

• 11 269 E— 02 
.138426-02 
•12488E-02 
.771416-03 
. 18 6 3 3 E -0 3 

— . 3898 8 E — 0 3 
-.893016-03 
-.117196-02 
-.1201 5E-02 
-. 10626E-02 
-.845206-03 
-.720346-03 
-.807346-03 
-.109446-02 
-. 15513E-02 

— • 21 572E-02 
-.283986-02 
- • 357 2 4 E — 02 
-.42746E-02 

— • 492 0 7E -0 2 
-.552086-02 

— .6041 5 E — 02 

— • 648 5QE-02 
-.686436-02 
- . 722 54E-02 
-.758296-02 

— • 79 37 5E— 02 
-.829046-02 

— • 86 433E— 02 
-.900086-02 
-.940236-02 
-.99074E-02 
- • 105 21 E-01 
-.112376-01 
- . 1203 9E-01 
-.129026-01 
- • 1378 IE-01 

— • 14481 E-01 
- • 148 60 E— 0 1 

— • 1 4900 E — 01 
-.145756-01 

— .13 8 64 E-0 1 
- . 127 50E-01 
-. 1146 7 E-0 1 
- . 104 96E— 0 1 

— • 10067 E-01 

— • 10358E-01 
-• 11434E — 01 
-.126086-01 

— • 1 308 0E— 0 1 

— • 12520E-01 
-.111616-01 
-.10285E-01 
-.104966-01 
-. 11228E-01 
-.111956-01 
- • 10 13 4 E-0 1 
-.847056-02 
-.66817E-02 

— • 543 4 9E-02 

— • 4905 1 E-02 
-.285076-02 


-.200006+01 

U 

-.607976-03 

• 9 75 94 E-0 3 
— . 5 4070F-02 

— • 75968 E— 02 
-.623856-02 
— .49950E-02 
-.549676-02 
-.784256-02 

— • 10522E-01 
-■ 11137E-01 

— .99 7 39 E— 02 
-.920356-02 
-.98645E-02 

— • 1 1 031 E-01 
-.1 1677E-01 
-.11255E-01 
-.993576-02 
-.810076-02 

— .60205 E-02 
- • 3 88 46E-02 
-.186416-02 
— . 1 0 896 E — 0 3 

.128436-02 

• 2 2 586E-02 
.278936-02 
.285206-02 
•25041E-02 
.186796-02 

• 1 0 666 E — 0 2 

• 2 2 3 31 E-03 
- . 5 38 60E-03 

— .1095 8 E— 0 2 
-.135726-02 
-.136196-02 
-•11820E-02 
-.889496-03 
- • 5 56 49E— 0 3 
-.254696-03 
-.559996-04 

.635896-05 
-.407606-04 
-. 16475E-03 
-.333036-03 
-.513006-03 
-.672066-03 
-.779 3 3E — 0 3 

— • 8 2 660 E —0 3 
-.821526-03 
-.771066-03 
-.678556-03 

— • 546 3 7E -0 3 
-.379946-03 
-.187646-03 

.187376-04 

.221396-03 

.393996-03 

• 5 02 1 7E-03 
. 56256E-03 

• 662 30E-03 
. 9 4022E-03 

• 1 59 4 4E— 0 2 
.259916-02 
.354936-02 
.400466-02 
.459616-02 
.637116-02 
•79775E-02 
•80820E-02 

• 9 1 1 9 4 E— 02 
.908206-02 
.961416-02 


— • 20000E + 0 1 

V 

— .10 9 89E— 02 
-.548536-02 
-.35 789 E — 02 

— .1091 9 E — 0 2 
- • 86 2 45 E — 0 3 
-.21811 E-02 
-.37291 E— 02 

— • 40 495 E-02 
-.244116-02 
-.84673E-03 
-.969 51 E— 0 3 
-.212376-02 
-.2574 3 E-02 
-.192906-02 
-.529626-03 

• 8 46 50 E-0 3 
.180596-02 
•23301E-02 

• 24 505E-0 2 
.217456-02 
.153466-02 
.620816-03 

-.478346-03 

— • 1 67 63 E-02 
-.295186-02 
-.407256-02 
-.489696-02 
-.538956-02 
-.55344E-02 
-.53328E-02 

— • 4 e040E-02 
-.403126-02 
-.324936-02 
- « 2 5505 E-02 
-.196846-02 
-.152456-02 
-.12297E-02 
-.108206-02 
-.102026-02 
— . 94078 E-03 
-.824646-03 

— • 66453E-03 
-.457296-03 
-.203626-03 

.89513 E-04 
•38231E-03 
.6443 7 E-0 3 
,870766-03 
.105946-02 
.121066-02 
.131526-02 
.136626-02 
. 13 590E-02 
.129026-02 
. 1 1 6 1 4 E-0 2 
.985706-03 

• 83 892 E-0 3 
•80447 E-03 
.924606-03 
.120176-02 
.139826-02 

• 1 2 302 E-02 
.784506-03 
.826906-03 
.145186-02 
.1161 0 E-0 2 

• 1 5079E-03 

• 40202 E-03 
-.671946-04 

.391986-03 



PHI IN VICINITY OF AIRFOIL 


J 

PHI ( 3.2) 

PHI <J» 2,2 ) 

PH I ( J > 1 > 2 ) 

PHI( J,l,l) 

PHI ( J» 2 j 1 ) 

PHI ( J 1 ) 

25 

-.50429F-01 

-•52140 E-01 

- • 5 1 8 0 1 E-0 1 

- • 5 1 80 1 E-01 

- • 4 752 6E-01 

-• 41340 E-01 

26 

-.56935E-01 

-•6158 4 F-01 

-•64478E-01 

- • 6 4478 F-01 

-• 5 561 3E-01 

- • 45505E-Q1 

27 

-•63086F-01 

-*72724 F-01 

- • 882 20E-01 

- • 9097 1 E-01 

- • 644 00 E-0 1 

-• 47826E-01 

28 

-.65524F-01 

-.75528E-01 

-.83473E-01 

-.67988F-01 

-• 54 5 56E-01 

-• 41802E-01 

29 

-.64117E-01 

-.72372E-01 

- • 8 1 5 47F-0 1 

- • 4 48 77E-01 

-• 3728 3E-01 

29235E-01 

30 

-.60009E-01 

—•661 59E-01 

- • 72 4 1 OE-0 1 

-.21256E-01 

- • 1736 3E-01 

- • 131 86E-01 

31 

- • 5 4407 E-01 

-•5 8722 E-01 

- • 628 7 8E-0 1 

• 12 8 49E-0 3 

• 1 4006E-02 

•29163E-02 

32 

- • 4 8 14 6 E-01 

- . 5 09 58 F-01 

- • 5 35 1 9F -0 1 

• 1 8 5 5 1 E-01 

•17959E-01 

• 17653E-01 

33 

-•41744 E-01 

-•43355F-01 

-•44717E-01 

• 34 070 E-0 1 

• 3234 5E-01 

• 30835E-01 

34 

-•35480E-01 

- • 3 6 1 0 4 E-0 1 

36506E-01 

• 4 7761 E-0 1 

• 4 5184E-01 

• 42741E-01 

35 

-•29518 E-01 

-•2 929 1 E-01 

-•28826E-01 

• 602 2 0E-01 

• 56829E-01 

• 53557E-01 

36 

- • 2 3992 F-01 

-•23072E-01 

-•21941E-01 

•71405E-01 

. 67303E-01 

• 63313 E-01 

37 

- • 1 8966 E-01 

- • 1 7464E-01 

- • 15779F-01 

• 8 1 3 30E-01 

•76612E-01 

•72012E-01 

30 

- • 14474E-01 

-.1246 5 F-01 

- • 1 02 69E -0 1 

• 8 994 5 E-01 

• 84747 E-01 

• 79659E-01 

39 

-•10 5 55 E-01 

-•8131 0 E-02 

-.5571 0E-02 

• 97274E-01 

• 91751E-01 

• 86284E-01 

40 

-•72124F-02 

-•43927E-02 

- • 141 59E-02 

. 10368 E + 00 

• 97776E-01 

• 91 949E-01 

41 

- • 44962 E-0 2 

-•13206E-02 

• 2 04 99E-02 

• 1090 3E +00 

•10278E+00 

• 96639E-01 

42 

-.24723E-02 

• 9 95 20 E-0 3 

• 470 6 9E-02 

• 1 1 3 2 6E + 00 

• 10670E+00 

• 10032E+00 

43 

- • 1 2069 E-02 

•24616E-02 

•64064F-02 

• 1162 7E + 00 

• 10948E + 00 

• 10294E + 00 

44 

-.75203F-03 

• 3 03 49 F-02 

•71366E-02 

• 1 1 796 E + 00 

• 11 102E + 00 

• 10444E+00 

45 

- • 1 1 747F-02 

• 2 6741 E-02 

• 709 1 2E-02 

• 1 1 81 5E + 00 

• 1 1 120E + 00 

• 10477E+00 

46 

- • 2620 IE-02 

• 1 01 67 E-02 

. 576 6 4E-02 

• 1 1 6 5 9E +00 

• 10984E + 00 

• 10388E + 00 

47 

-.52134F-02 

- • 2 62 96 E-02 

• 626 6 8E-03 

• 11231E+00 

• 10672 E + 00 

• 101 82E + 00 

48 

- • 8 4 36 3 E-02 

-•74551 E-02 

- • 727 5 9E-02 

•10440E+00 

• 10214E + 00 

•98992E-01 

49 

-• 11408 E-01 

-•11 2 55 E-01 

-•11746E-01 

• 99933 E-01 

• 98 504 E-01 

• 96324E-01 


LaJ 



I 

JMX 

KMX 

ISUP 

NS 

0 PHI MX 

DA 

U F AC 

EU 

1 

71 

17 

0 

4 

-.544586-01 

0. 

. 10000E+C1 

. 10000E+01 

2 

71 

18 

0 

4 

-• 30782E-01 

0. 

• 10000E+0 1 

. 10000E+01 

3 

71 

17 

0 

4 

-.191 30E-01 

0. 

• 10000E + 0 1 

• 1 OOOOE+Ol 

4 

68 

3 

0 

4 

-.13160E-01 

0. 

• 10000E +0 1 

• lOOOOE+Ol 

5 

6e 

2 

0 

4 

-• 96538E— 02 

0. 

. lOOOOE+Ol 

• 10000E + 01 

6 

67 

7 

0 

4 

-.769716-02 

0. 

• 10000E+01 

• lOOOOE+Ol 

7 

66 

7 

0 

4 

-.639596-02 

0. 

• 10000E+01 

• lOOOOE+Ol 

8 

65 

5 

0 

4 

-•52869E— 02 

0. 

• 1 OOOOE+Ol 

•lOOOOE+Ol 

q 

6 4 

2 

0 

4 

-.444686-02 

0. 

• 10000E +0 1 

*10000 E+ 01 

10 

63 

2 

0 

4 

— . 37217E— OH 

0. 

. 10000E+01 

•lOOOOE+Ol 

li 

62 

2 

0 

4 

- • 3 101 5E— 02 

0. 

* 10000E+01 

.100006+01 

12 

61 

2 

0 

4 

-.25836E-02 

0. 

• 10000E+01 

« 10Q00E+01 

13 

60 

2 

0 

4 

-.2 15 78 E— 02 

0. 

.lOOOOE+Ol 

• 10000E + 01 

14 

59 

2 

0 

4 

-.130956-02 

0. 

•10000E+01 

*1 OOOOE+Ol 

15 

59 

-6 

0 

5 

-.15884E-02 

0. 

. 10000E+01 

•1 OOOOE+Ol 

16 

58 

-6 

0 

4 

-.140436-02 

0. 

. 10000E + 01 

* 1000QE+01 

17 

57 

-6 

0 

5 

-.123936-02 

0. 

.100006+01 

• lOCOOE+Ol 

18 

56 

-6 

0 

4 

-.10929E-02 

0. 

. 10000E + 01 

* lOOOOE+Ol 

10 

55 

-6 

0 

5 

— . 96379E-03 

0. 

• 10000E+01 

• 1 OOOOE + Ol 

20 

54 

4 

0 

5 

-.85178E-03 

-* 38431E-03 

. 10000E+01 

• 1 OOOOE+Ol 

21 

53 

6 

0 

5 

- • 7 5623 E— 03 

-.74701E-03 

•10000E+01 

•1 OOOOE+Ol 

2? 

52 

6 

0 

4 

—.672736—03 

-.107106-02 

•10000E+01 

•10C00E+01 

23 

51 

6 

0 

4 

— .60236E— 03 

-♦ 1 3531E-02 

. 10000E+01 

• 1 OOOOE+Ol 

24 

50 

6 

0 

4 

- • 5461 4E— 03 

-.15931E-02 

» 10000F+01 

• 1 OOOOE + Ol 

25 

49 

6 

0 

4 

- • 50069E— 03 

-.179336-02 

. 1000CE+01 

•lOOOOE+Ol 

26 

48 

6 

0 

4 

— .4 6169E-03 

-.19570E— 02 

. 10000E+01 

. lOOOOE+Ol 

27 

48 

8 

0 

4 

-• 42773E-03 

— . 2 0877E-02 

. 10000E+01 

• 1 0000E+ 01 

28 

47 

8 

0 

4 

— • 40090E-03 

- ■ 2 1895E-02 

. 10000E + 01 

•1G000E+01 

20 

46 

8 

0 

4 

-.374 5 7 E— 03 

— • 22666E-02 

.100006+01 

• 1 OOOOE+Ol 

3 C 

45 

8 

0 

4 

-.34935E-03 

-• 23230E-02 

• lOGOOE+Ol 

• 1 OOOOE + Ol 

31 

4 4 

8 

0 

4 

- • 3 24 8 3E-03 

- • 2 3624E-Q 2 

•10000E+01 

•lOOOOE+Ol 

32 

43 

8 

0 

4 

-« 30068E-03 

-.23878E-02 

. 10000E + 01 

• 10000E + 01 

33 

42 

8 

0 

4 

-.277566-03 

- • 2 4024E-02 

. 10000E+G1 

.100006+01 

34 

4 3 

1 

0 

4 

- • 2 5624E— 03 

-• 24O08E-O2 

• 100 00 E + Gl 

•1C000E+01 

35 

42 

1 

0 

4 

- • 2 3976E — 03 

-.24092E-02 

• 1 0000E + 01 

.100008+01 

36 

41 

1 

0 

4 

-.22366E— 03 

-.24054E-02 

. 10000E + 0 1 

• 10000F+01 

37 

4 0 

1 

0 

4 

-.20834E-03 

-.23990E-02 

. 10000E+0 1 

• 1 OOOOE+Ol 

3 P 

39 

1 

0 

4 

— • 19 40 3E— 0 3 

-.2 3909E-02 

•10000E+01 

• 1 OOOOE+Ol 

30 

38 

1 

0 

4 

— .18 090E — 03 

— • 2 3621E — 02 

•10000E+01 

. lOOOOE+Ol 

<♦0 

38 

1 

0 

4 

-.16990E-03 

-.237316-02 

. 1 OOOOE+Ol 

• lOOOOE+Ol 

41 

37 

1 

0 

4 

— .15988E — 03 

-.23644E-02 

• 1 OOOOE + O 1 

• 1 OOGOE + 01 

42 

36 

1 

0 

4 

-.15090E-03 

— . 2 3562E-02 

• 10000E +01 

• 1 OOOOE + Ol 

43 

35 

1 

0 

4 

— .1424 IE— 03 

-.234876-02 

• 100 OOE +0 1 

• lOOOOE+Ol 

44 

34 

1 

0 

4 

- • 1 3 307E-03 

-.2 3420E— 02 

• 1 OOOOE + Ol 

• 1000CE+01 

45 

32 

1 

0 

4 

— . 12390E— 03 

-.233636-02 

. lGOQQE+Ol 

• 1 0000 E + 01 

46 

32 

1 

0 

4 

-.10928E-03 

— • 23316E-02 

. 10000E+01 

•1 OOOOE+Ol 

47 

43 

-1 

0 

4 

-* 97027E — 04 

-.23282E-02 

•10000E+01 

. 1 OOOOE + Ol 

4 8 

48 

-1 

0 

3 

-.17427E-03 

- • 2 293 5E-02 

•99890E+00 

• 5 5 66 2 E-Ql 

49 

47 

-1 

0 

3 

-• 12114E— 03 

- * 22 61 6E -0 2 

.99890E+00 

.55662E-01 

5C 

46 

-1 

0 

3 

-• 14912E— 03 

- . 2 2527E-02 

•99890E+00 

• 55662E-01 

51 

45 

-2 

0 

3 

-.10182E— 03 

— • 22417E-Q2 

• 998 90 E +00 

• 55662E-01 

52 

48 

-1 

0 

3 

-.12091 E -03 

-.22425E-02 

• 99890E + 00 

. 55662E-01 

53 

43 

-2 

0 

3 

-.89928 E— 04 

-.22400E-02 

. 99890F+0C 

• 5 5662 E-01 

54 

48 

-1 

0 

3 

-.2010 IE-03 

-.22087E-02 

• 997 69 E +00 

•55611E-01 

55 

47 

-1 

0 

3 

-.117496-03 

-.21766E-02 

• 997 69F + 00 

.556116-01 

56 

48 

-1 

0 

3 

-.16395E-03 

— .21727E-02 

•99769E+00 

•55611E-01 

57 

45 

-2 

0 

3 

- . 9 49 7 5E— 04 

-.21636E-02 

•99769E+00 

.556116-01 

58 

48 

-1 

0 

3 

-* 12745E-03 

-.2 1690E-02 

•99769 E+OC 

.556116-01 

59 

43 

-2 

0 

3 

- . 82 32 7E— 04 

-.216876-02 

• 99769 E + 00 

.556116-01 

60 

5 C 

1 

0 

3 

-.7617 5E-04 

-.21 949E— 02 

• 998 30E + 00 

• 56 199E-01 

61 

41 

-2 

0 

3 

-.602346-04 

— .22 157E-02 

.998306+00 

.561996-01 

62 

40 

-2 

0 

3 

-.64215E-04 

-.22293E-02 

.99830E+00 

.561996-01 

63 

46 

-1 

0 

3 

- * 63 53 2 E-04 

- • 2 2 41 2E — 02 

•99830E+00 

. 56199E-01 

64 

36 

-2 

0 

3 

-.59048E-04 

-.22498E-02 

• 99830E+00 

.561996-01 


STEPS 

65 - 

195 

omitted 






32 



196 

49 

-1 

0 

4 

- • 2 04 53E-04 

-.255516-02 

.10014F+01 

. 54774E — 01 

197 

50 

-1 

0 

4 

. 14432E-04 

-.255206-02 

.10014E+01 

.547746-01 

198 

50 

-1 

0 

4 

- • 1 354 5E— 04 

- • 2 5506E— 02 

* 1 00 1 4 E +01 

.547746-01 

199 

20 

4 

0 

4 

— . 96440E-05 

-.254766-02 

• 10014E + 01 

• 54 77 4E-01 

2 CO 

46 

-1 

0 

4 

.72986E-04 

-.25744E-02 

• 10024E + 01 

• 54666E-01 

201 

48 

-1 

0 

4 

— •6618 8 E— 04 

-.25988E-02 

• 10024E+01 

• 5 4666E — 0 1 

202 

48 

-1 

0 

4 

• 58745E— 04 

-.260496-02 

• 1 00 24E +0 1 

• 54 666 E-0 1 

203 

48 

-1 

0 

4 

-.485086-04 

-.261306-02 

. 10024E+01 

• 5 466 6E — 01 

204 

46 

-1 

0 

4 

• 4 0 1 3 5E— 04 

-.261226-02 

. 10024E + 0 1 

. 54666E-01 

205 

40 

-1 

0 

4 

— . 32088 E— 04 

-• 26144E— 02 

• 10024E + 01 

.546666-01 

206 

49 

-1 

0 

4 

• 2 591 9E-04 

— .26113E-02 

.10024E+01 

.54666E-r01 

207 

49 

-1 

0 

4 

— • 2091 5E— 04 

-.261076-02 

• 1 00 2 4E + 0 1 

• 54 666E — 0 1 

208 

50 

-1 

0 

4 

• 1 64 37 E-04 

-.260716-02 

•10024E+01 

• 5 4 66 6E — 0 1 

209 

50 

-1 

0 

4 

— • 1 3 61 7E — 04 

-.260536-02 

•10024E+01 

. 54666E-01 

210 

51 

-1 

0 

4 

• 1 032 IE — 04 

— .26018E-02 

• 10024E + 0 1 

• 5 4 66 6E — 01 

211 

19 

4 

0 

4 

— • 91 039E-05 

- • 2 5996E— 02 

.10024E+01 

• 5 4666 E—0 1 

212 

4 P 

-1 

0 

4 

— • 34696E— 04 

- • 2 5823E— 02 

• 100 19E + 0 1 

• 54700E — 01 

213 

48 

-1 

0 

4 

.298556-04 

-.256686-02 

.100196+01 

•54700E-01 

214 

46 

-1 

0 

4 

—.2861 7E— 04 

-• 25603E— 02 

.100196+01 

• 5 4700E — 01 

215 

4 8 

-1 

0 

4 

• 21 697E-04 

-.255336-02 

• 1 00 19E + 0 1 

• 54700E— 01 

216 

48 

-1 

0 

4 

- • 2 0077E-04 

-. 2 5510E— 02 

•10Q19E+01 

• 54700E— 01 

217 

48 

-1 

0 

4 

• 1 3892E— 04 

— • 2 5475 E— 02 

.100196+01 

.547006-01 

218 

49 

-1 

0 

4 

-. 135 3 OE -04 

-.254696-02 

.100196+01 

.547006-01 

219 

49 

-1 

0 

4 

.8514BE-05 

-.254536-02 

•10019E+01 

• 54700E— 01 

220 

29 

2 

0 

4 

— • 16609E — 04 

— .2 5 526E-02 

• 1 00 2 1 E +01 

• 54 7 52 E — 01 

221 

48 

-1 

0 

4 

— • 1 3042E-04 

-.2 5 588E — 02 

•10021E+01 

• 5 4752 E — 01 

222 

4 e 

-1 

0 

4 

• 1 022 7 E— 04 

-.25612E-02 

•10021E+01 

• 5 47 5 2 E-0 1 

223 

46 

-1 

0 

4 

- • 1 02 10E— 04 

— • 2 5638E— 02 

• 10021E + 01 

.547526-01 

224 

48 

-1 

0 

4 

. 6 9 0 3 5 E —0 5 

- • 2 5 646E — 0 2 

.100216+01 

• 5 475 2E-01 

225 

29 

2 

0 

4 

— .2018 3 E— 04 

-.2 5736E-02 

.100246+01 

• 5 47 3 1 E — 01 

226 

4 6 

-1 

0 

4 

— • 1540RE — 04 

— • 2 5810E-02 

• 1 00 2 4E +0 1 

. 54731E-01 

227 

46 

-1 

0 

4 

• 1 31 3 2E— 04 

— .2 5836E-02 

.10024E+01 

• 5 4 7 31 E — 01 

228 

48 

-1 

0 

4 

-.11849E-04 

— • 2 5862E— 02 

. 10024E+01 

.547316-01 

229 

48 

-1 

0 

4 

• 90269E-05 

- . 2 5864E — 02 

. 10024E+01 

•54731E-01 

230 

29 

2 

0 

4 

— • 22689E — 04 

-.259576-02 

.10027E+01 

.547296-01 

231 

48 

-1 

0 

4 

— • 15P39E — 04 

- • 2 603 IE-02 

.100276+01 

. 54 729 E — 01 

232 

46 

-1 

0 

4 

• 1433 5E — 04 

— .2 605 3E — 02 

. 10027E+01 

.547296-01 

233 

48 

-1 

0 

4 

—.121 96E— 04 

- • 2 6074E— 02 

•10027E+01 

• 5 4 7 2 9 E— 01 

234 

48 

•1 

0 

4 

• 1 0012E — 04 

- . 2 6070E— 02 

•10027E+01 

.547296-01 

235 

4« 

-1 

0 

4 

— • 82390E-05 

— • 2 6072E — 02 

. 10027E+01 

.547296-01 

236 

50 

-2 

0 

4 

- • 71 206 E—0 5 

— .2 601 6 E— 02 

• 10025E + 01 

.547588-01 

237 

22 

2 

0 

4 

— • 57290E— 05 

-.25969E-02 

• 10025E+01 

• 54758E-01 

2 3 P 

48 

-1 

0 

4 

— • 5 996 5E — 0 5 

- • 2 5 94 2 E— 02 

•1Q025E+01 

. 54758E-01 

239 

2 C 

2 

0 

4 

-•46322E-05 

-.25917E-02 

• 1 0025E +01 

•54758E-01 

240 

48 

-1 

0 

4 

44780E— 05 

-.2 5901E-02 

• 100 2 5E +0 1 

•54758E-01 

241 

46 

-1 

0 

4 

• 0962 1 E — 05 

— . 2 5907E — 02 

.1Q026E+01 

• 5 47 8 1 E — 0 1 

242 

48 

-1 

0 

4 

- • 8 3 757E-0 5 

-.259186-02 

• 1 00 2 6 E +0 1 

•54781E-01 

243 

48 

-1 

0 

4 

• 6 4 344 E —0 5 

— .25914E-02 

. 10026E+01 

• 54781E-01 

244 

48 

-1 

0 

4 

— . 58743E — 05 

-.25915E-02 

• 10026E+01 

.54781 E -01 

245 

46 

-1 

0 

4 

• 4109 3E— 05 

— • 2 591 0 E-02 

•10026E+01 

.54781E-01 

246 

48 

-1 

0 

4 

70167E-05 

— .25898 E — 02 

. 10026E + 01 

.547926-01 

247 

4 8 

-1 

0 

4 

.51 50 2 E— 05 

-.25682E-02 

. 1 00 2 6 E + 0 1 

. 5479 2 E — 01 

24 R 

48 

-1 

0 

4 

-.5091 1 E — 0 5 

-.25 878E-02 

•10026E+01 

. 5 479 2E-01 

249 

48 

-1 

0 

4 

. 33549E-05 

-. 25869E-02 

• 10026E +0 1 

• 54 79 2 E-01 

250 

48 

-1 

0 

4 

-.35104E-05 

— . 2 5 867E-02 

• 10026E +01 

• 54 79 2 E — 01 

251 

48 

-1 

0 

4 

. 3574 IE— 05 

- • 2 58 67E-02 

.100266+01 

.547976-01 

252 

48 

-1 

0 

4 

— • 37409E— 0 5 

-.25870E-02 

.100266+01 

• 54797E-01 

253 

48 

-1 

0 

4 

• 23641E — 05 

-.25B60E-O2 

. 1 0026E +01 

• 5 4 79 7 E — 01 

254 

48 

-1 

0 

4 

-.264096-05 

— .2 5 869E-02 

• 10026E+01 

• 54 797E-01 

255 

46 

-1 

0 

4 

.141196-05 

— • 25867E— 02 

.100266+01 

• 5 479 7E-01 

256 

48 

-1 

0 

4 

- • 2 6 13 2E — 05 

-.258646-02 

• 100 26 E + 0 1 

• 54 799E — 0 1 

257 

48 

-1 

0 

4 

■ 1 4239E— 05 

— • 2 5860E — 02 

.100266+01 

.547996-01 

258 

48 

-1 

0 

4 

-.19326E-05 

— .25859E— 02 

•10026E+01 

• 5 4799E — 01 

259 

21 

7 

0 

4 

-• 93532E — 06 

— . 2 5856E-02 

• 10026E+0 1 

• 5 4799 E — 01 

260 

48 

-1 

0 

4 

— • 1 41Q9E— 05 

— .2 5855 E— 02 

• 10026E+01 

. 54799E-01 
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u> 

4 > 


EQUIVALENT FREF AIR (U AND V APE PFR TURB AT I QMS FROM U INF BUT NORMALIZED BY UREF) 


AMINF = .68906, UINF/URFF = 1.00258> CPStAR * -.82530F+00, OFAC = 1.00393 ALPHA CORRECTION * -.25855E-02 


X (J-.5) 

■ YMOD 

CP ( INF ) 

U 

J 

X ( J ) 

V 

DELTA PHI 

DELTA V 

-.02295 

0.00000 

.89500 

-.53133 

27 

.00200 

.79856 




0.00000 

.80900 

-.48104 



-.47098 

-.00247 

1.27754 

.02695 

.03975 

-.61591 

.46041 

28 

.05190 

.25735 




-.02384 

.14563 

-.01006 



-.27652 

.02105 

.53387 

.07685 

.05255 

-.92393 

.46095 

29 

.10180 

.14646 




-.03760 

-.21455 

.13314 



-.17905 

.03740 

.32551 

.12675 

.05984 

-.98364 

.46030 

30 

.15170 

.07574 




-.04651 

-.34113 

.17753 



-.11720 

.05191 

.19294 

.17665 

.06361 

-.89566 

.42122 

31 

.20160 

.02900 




-.05235 

-.37360 

.18552 



-.07571 

.06367 

.10471 

.22655 

.06506 

-.76973 

.36207 

32 

.25150 

-.01054 




-.05612 

-.37343 

.18209 



-.04423 

.07266 

.03369 

• 2 764 5 

.06453 

-.64438 

.30522 

33 

.30140 

-.03296 




-.05832 

-.35205 

.17098 



-.02128 

.07935 

-.01167 

, 32635 

.06289 

-.56172 

.26813 

34 

.35130 

-.04877 




-.05930 

-.32874 

.15913 



-.00319 

.08479 

-.04558 

.37625 

.06047 

-.50541 

.24318 

35 

.40120 

-.06475 




-.05953 

-.30620 

.14845 



.01431 

.08952 

-.07907 

.42615 

.05724 

-.44742 

.21751 

36 

.45110 

-.07868 




-.05882 

-.27201 

.13251 



.02683 

.09376 

-.10551 

.47605 

.05333 

-.39036 

.19214 

37 

.50100 

-.09096 




-.05749 

-.24073 

.11799 



.03718 

.09746 

-.12814 

.52595 

.04880 

-.33194 

.16583 

38 

.55090 

-.10046 




-.05564 

-.21236 

.10486 



.04760 

.10050 

-.14814 

.57585 

.04380 

-.27570 

.14001 

39 

.60080 

-.10562 




-.05326 

-.17752 

.08853 



.05354 

.10307 

-.15915 

.62575 

.03854 

-.23531 

.12150 

40 

.65070 

-.11374 




-.05060 

-.15397 

.07757 



.06220 

.10527 

-.17593 

.67565 

.03298 

-.10908 

.10018 

41 

.70060 

-.12089 




-.04750 

-.12426 

.06367 



.07036 

.10709 

-.19125 

.72555 

,02686 

-.14094 

.07764 

42 

.75050 

-.12734 




-.04400 

-.08984 

.04739 



.07775 

.10660 

-.20509 

.77545 

.02053 

-.00983 

.05330 

43 

.80040 

-.13265 




-.04013 

-.04989 

.02813 



.08246 

.10985 

-.21510 

.82535 

.01392 

-.03554 

.02691 

44 

.85030 

-.13697 




-.03603 

-.01023 

.00861 



.08459 

.11077 

-.22156 

.87525 

.00711 

.02491 

-.00319 

45 

.90020 

-.13911 




-.03182 

.02181 

-.00705 



.09361 

.11096 

-.23272 

.92515 

.00018 

.09331 

-.03803 

46 

.95010 

-.14224 




-.02716 

.07518 

-.03280 



.11646 

.11070 

-.25870 

.97505 

-.00690 

.19517 

-.09254 

47 

1.00000 

-.12821 




-.02136 

.21924 

-.10918 



.08352 

.11153 

-.21173 


-.0132P 

-.01720 


1.02495 



EXPERIMENTAL CP BASED ON CORRECTED AMINF, CL 


222A0 


X CP(INF) 


.02695 
.07685 
.12675 
.17665 
.22655 
. 2 76 A 5 
. 32635 
.37625 
. A2615 
. A 760 5 
.52595 
. 57585 
.62575 
.67565 
.72555 
. 775 A 5 
.82535 
.67525 
.92515 
.97505 

RMS DELTA CP = 


— . 61 2 A8 
. 1AAA2 

— » 92 A 30 

-.21909 
-.98981 
-• 3A629 
-.90680 
-.37913 
-.77983 
-.37918 
-.65113 
-.35871 
-.56927 
-. 33A71 
-. 513A0 
-.31232 
-. A5559 
-.27822 
-.39861 
— . 2 A70 6 
-. 3A020 
-. 2188A 
-.28396 
-. 18A09 

— . 2 A 3 70 

-.16073 
-. 197A5 
-.13117 
-. 1A925 
-.09689 
-.09803 
-.05705 
- , 0A3 57 
-.01750 
.01712 
. 01A3 5 
. 0858A 
.06759 
.16829 
.21221 

. 70930E-02 
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UPPFR 

SURFACE 

LOWER 

SURFACE 

X 

CP(INF) 

X 

CP(INF) 

.00022 

1.08939 

.00022 

1.08939 

.01436 

-.32484 

.01172 

.50168 

.02822 

-.61813 

.02434 

.18546 

*06297 

-.83806 

.04994 

-.09600 

.07798 

-.92659 

.07417 

-.20924 

.10296 

-.98393 

.09855 

-.29178 

.15231 

-.95264 

.14860 

-.36580 

.20155 

-.85013 

.19846 

-.38267 

.25228 

-.70696 

.24811 

-.37188 

.30200 

-.60486 

.29807 

-.34786 

.35197 

-.53853 

.34861 

-.32497 

.40145 

-.48644 

. 39795 

-.29966 

.45146 

-.42429 

.44826 

-.26199 

.50149 

-.37157 

.49806 

-.23617 

. 55082 

-.30872 

.54869 

-.20280 

.60089 

-.26364 

.59798 

-.17166 

.65075 

-.22182 

.64853 

-.15034 

.70067 

-.17255 

.69885 

-.11356 

.74982 

-.12580 

.74836 

-.08144 

.80091 

-.06963 

.79659 

-.03607 

.85039 

-.01527 

.84924 

-.00326 

.89945 

.05028 

.89914 

.03500 

.94876 

.12668 

.94766 

.11670 

.96750 

.22079 

.98750 

.26259 


I 

J MX 

KMX 

ISUP 

NS 

DPHIMX 


10 

65 

-1 

0 

4 

• 73865E-04 


2 0 

4 6 

18 

0 

4 

•16297E-04 


30 

37 

19 

0 

4 

• 70191E-05 


40 

22 

-3 

0 

4 

-. 40991E-05 


50 

15 

-3 

0 

4 

-.17667E-05 


55 

11 

-4 

0 

4 

-.84840E— 06 

0 . 
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MODEL PFRTIIRBATION (U AND V ARE PERTURBATIONS FROM UINF 8UT NORMALIZED BY UREFI 


MJ-.5) 

YMOD 

CP< INF) 

U 

J 

-•02295 

0.00000 

.89752 

-.53394 

27 


0.00000 

.80539 

-.47800 


.02695 

.03982 

-.60690 

.45661 

20 


-.02376 

.14161 

-.00907 


•07605 

.05270 

-.91793 

.45053 

29 


-.03746 

-.21731 

.13419 


.12675 

.06005 

-.98236 

.46790 

30 


-.04631 

-.34173 

.17763 


•17665 

.06387 

-.89860 

.42262 

31 


-.05210 

-.37290 

.18510 


.22655 

.06534 

-.77237 

.36328 

32 


-.05583 

-.37254 

.18164 


.27645 

.06484 

-.64479 

.30539 

33 


-.05802 

-.35247 

.17079 


.32635 

.06323 

-.56277 

.26859 

34 


-.05905 

-.32801 

.15878 


.37625 

.06082 

-.50674 

.24377 

35 


-.05919 

-.30521 

.14799 


.42615 

.05762 

-.44889 

.21816 

36 


-.05845 

-.27087 

.13198 


.47605 

.05372 

-.39187 

.19281 

37 


-.05710 

-.23951 

.11741 


.52595 

.04921 

-.33341 

.16650 

38 


-.05524 

-.21112 

.10428 


• 575B 5 

.04422 

-.27713 

.14067 

39 


-.05286 

-.17628 

.08794 


.62575 

.03897 

-.23672 

.12216 

40 


-.05019 

-.15273 

.07697 


.67565 

.03331 

-.19042 

.10002 

41 


-.04709 

-.12305 

.06309 


.72555 

.02730 

-.14219 

.07826 

42 


-.04359 

-.06071 

.04683 


.77545 

.02096 

-.09097 

.05387 

43 


-.03972 

-.04085 

.02760 


.82535 

.01435 

-.03653 

.02743 

44 


-.03562 

-.00933 

.00815 


.87525 

.00753 

.02414 

-.00277 

45 


-.03142 

.02252 

-.00744 


.92515 

.00059 

.09275 

-.03771 

46 


-.02677 

.07568 

-.03300 


.97505 

-.00650 

.19475 

-.09228 

47 


-.02099 

.21958 

-.10939 


1.02495 

-.01289 





-.01685 


X(J) 

V 

DELTA PHI 

DELTA V 

.00200 

. 80010 
-.47746 

-.00275 

1.27756 

.05190 

.25864 

-.27525 

.02049 

.53389 

.10180 

.14775 

-.17777 

.03667 

.32552 

.15170 

.07673 

-.11622 

.05115 

.19296 

.20160 

.02964 

-.07509 

.06301 

.10473 

.25150 

-.01012 

-.04383 

.07207 

.03371 

.30140 

-.03244 

-.02079 

.07879 

-.01165 

.35130 

-.04829 

-.00274 

.08427 

-.04556 

.40120 

-.06434 

.01470 

.08905 

-.07904 

.45110 

-.07834 

.02715 

.09335 

-.10549 

.50100 

-.09068 

.03743 

.09711 

-.12811 

.55090 

-.10024 

.04787 

.10021 

-.14611 

.60080 

-.10546 

.05367 

.10284 

-.15912 

.65070 

-.11364 

.06226 

.10510 

-.17590 

.70060 

-.12085 

.07037 

.10698 

-.19122 

.75050 

-.12736 

.07769 

.10855 

-.20505 

.80040 

-.13274 

.08233 

•10936 

-.21507 

.85030 

-.13713 

.08440 

.11082 

-.22153 

.90020 

-.13932 

.09337 

.11106 

-.23269 

.95010 

-.14246 

.11620 

.11003 

-.25067 

1.00000 

-.12854 

.08315 

.11168 

-.21169 



UJ 

00 


VfinCITIFS 0M CENTFR LIMF. 


XV 

-• 39083F+01 
-.34453E+01 
-• 3 0 5 34E + 0 1 
2 71 76E + 0 1 
24265E+01 
-.2171RE+01 
19470F+01 
-• 17473E+01 
1 5605E+O1 
- • 1 4077E +01 
1 262 1 E+01 
-• 1129PF+01 
“• 10090E+01 
89 R? BE + 00 
-. 79640E+00 
-.70236E+00 
-•61528E+00 
53443F+00 
45914E+00 
38B0BF+OO 
-• 32315E+00 
26153E+00 
-• 20365E+00 
-• 14917E+00 
-•97800 F -01 
-. 47900E-01 
•20000E-02 

•51900E-01 

. 1 0 18 CE + 00 

. 15170F+00 

• 2 01 60 E +00 
. 25150E+00 
• 30140E+00 
. 35130E+00 

• 40 1 2 OE + 00 


VF 

•30985E-02 

• 3 4901 F -02 
•39067F-02 
•43519E-02 

• 4 82 88E-02 
• 53414E-02 
. 58940E-02 
•64916F-02 
.71402F-02 
•78470F-02 
• 86206F -02 
•94715F-02 
. 10412F-01 
•11459E-01 
. 12632E-01 
•13950F-O1 
•15469E-01 
.172116—01 
•19246F-01 

• 21 661E-01 
•24504F-O1 

• 2 8 207E -01 

• 32 830F -01 

• 38 98 5 F -01 
•47528E-01 
•60424E-01 

• 800 10E +00 
-.4774 6E +00 

• 2 5 8 64 E +00 
-•27525E+00 

• 14 775E +00 

17777E+00 
• 76734E-01 
-•11622E+00 

• 296 44E-01 
-.75085E-01 
10120E-01 
-.43830E-01 

-• 324 40E—01 
-.20790F-01 
-• 48292E-01 
- . 2 735 6E-02 
- • 64 3 36E -01 
•14703E-01 
- • 78 336E-01 
• 2 7149E-01 


VW 

- • 42 382F-02 
-.41572E-02 
-.41839E-02 
-.42075E-02 
-•42212E-02 
-.42316E-02 

- • 42 379E-02 
-.4236RE-02 
-•42262E-02 
-.42066F-02 
-•41796E-02 
-•41 472E-02 
-•41111E-02 
- • 40 730E-02 

- • 4034 IE-02 

- • 39955E-02 
- • 39 577E-0 2 
-• 39215E-02 

- • 38 872F-02 
-• 30555E-O2 

- • 3 0269E-O2 

- • 38023 E-02 
- • 37 832F-0 2 
- • 37719E-02 

- • 3773 1 E-02 
- • 37 96 8E-0 2 
-•41126E-02 

-•41 146E-02 

- • 38705E-02 

-• 38701E-02 
-• 38704E-02 
-• 38681E-02 
-• 35711E-02 
-• 35667E-02 
-.32258E-02 
-.32197E-02 

- • 30021 E-02 

-• 29942E-02 
- « 30957E-02 
-• 30863E-02 

- • 30545E-02 

- • 30437E -02 
- . 2 9940E-02 
-.298 18E-02 

- • 292 5 7E-0 2 

- • 2 9 12 2 E -02 


4 5 1 1 OF + 00 


XU 


UF 


uw 


-•36772E+01 
-.32497E+01 
-• 28857E+01 
-.25722E+01 
-.22 993E+01 
- • 20595E +01 
-* 18472E + 01 
16580E + 01 
-.14882E+01 
-•1334 9E +01 
- • 11960E +01 
-•10694E+01 
- • 953 67E + 00 
- • 847 36F + 00 
- • 749 40E +00 
-.65884E+00 
-• 57487E + 00 

- • 49680E + 00 
-• 42402E+00 
- • 35 603 E + 00 

- • 2923 5E +00 
- • 2 3260E + 00 
- • 1764 IE + 00 

12349E + 00 
- • 72850E-01 

- • 22 9 50E-01 

•26950E-01 
• 76850E-01 
•12675E+00 

• 1766 5E + 00 

• 226 5 5E +00 

• 27645E + 00 
•32635E+00 

• 3762 5E +00 

• 42615F + 00 


-• 16325E-02 
-.19479E-02 
-.23359E-02 
- • 2 7996E-02 
-• 3 3465F-02 
-• 39870E-O2 
-.47383F-02 
-.56168E-02 
-.66470E-02 
-.78509E-O2 
- • 92900 E-02 
-•10988E-01 
-.13016E-01 
- • 1 54 5 3E-0 1 
- • 1 8409E-01 
- • 2 20 30E-01 
-.26521E-01 
“• 321 78E-01 
-.39442E-01 
- • 4 9000E-0 1 
- • 6 2000E-01 
-•80532E-01 
-.10894E+00 
-•15790E+00 
-.25723E+00 
-.53394E+00 
- • 4708 OE + 00 

• 4 5 661 E + 00 

— • 90651 E-02 
•45853E+00 

• 1 34 1 9E +00 

• 46790E +00 

• 1 7763E + 00 

• 4 22 62 E +00 

•18510E+00 

•36328E+00 

•18164E+00 

•30539E+00 

. 17079E+00 
•26859E+00 

•15878E+00 

•24377E+00 

• 1 4799E + 00 

• 2 1 8 1 6F + 00 

• 13198E + 00 


• 20 171E-05 
- • 16033 E-03 
-.34861 E-03 
- • 4 7 758E-0 3 

- • 52281 E-03 
-•49086E-03 

- .393416-0 3 
-.24121 E-03 

- • 44745E-04 

• 18 5 38 E-03 

• 43 889E-03 

• 70 64 3 E-03 

• 97 995E-03 

• 1 2 528 E-02 

• 15 197E-02 

• 1 7767E-02 

• 202 10E-02 

• 22 505E-02 
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Wall pressure range 



Figure 1.- Schematic representation of upper half of computational grid. Lower 
half is mirror image about x axis. 
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- 1.2 -.8 -.4 0 .4 .8 1.2 1.6 2.0 

X/C 

NflCR 0012, M=.701, RLPHR=2 .025 , RN=3 .Q46E+06 . 

SIGU= .00088 SIGV= .00052 SIGR= -00102 
RRIRTTM 83/07/22. 12.27.10. 


Figure 2.- First frame of plotted output for sample case. 


42 



x/c 

NflCfl 0012, M=.701, flLPHfl=2 .025 , RN:=3 -046E+06 
MINF= .68906 CL= .22248 SIGCP= .00709 
RRIRTTM 83/07/22. 12-27.11. 


Figure 3.- Second frame of plotted output for sample case. 
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NflCfl 0012. Mr. 701. flLPHRr2.025. RNr3 -046E+06 . 
SIGUr .00557 SIGVr .00181 SIGRr .00586 
RRIPTTM 83/07/22. 12-23.31. 


Figure 4. First frame of plotted output from unified 
procedure using sample case data. 


x/c 

NRCR 0012, Mr. 701, RLPHRr2 .025 , RNr3.046E+06 
MINFr .69757 CLr .22189 SIGCPr .01908 
RRIRTTM 83/07/22. 12.23.33. 


Figure 5.- Second frame of plotted output from unified 
procedure using sample case data. 
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